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Abstract 


The  purpose  of  this  research  was  to  model  and  simulate  the 
performance  of  a  digital  phase  shift  keyed  satellite  modem.  The 
probability  of  bit  error  (Pt)  at  different  levels  of  energy  per 
bit  to  noise  power  ratio  (Et/No)  was  the  performance  measure.  The 
channel  was  assumed  to  contribute  only  additi^’^e  white  Gaussian 
noise. 

A  second  order  Costas  loop  performs  demodulation  in  the 
modem  and  was  the  key  part  of  the  simulation.  The  Costas  loop 
with  second  order  Butterworth  arm  filters  was  tested  by  finding 
the  response  to  a  phase  or  frequency  step.  The  Costas  loop 
response  was  found  to  be  in  agreement  with  theoretical 
predictions  in  the  absence  of  noise.  Finally,  the  effect  on  Pt 
of  a  rate  1/2  constraint  length  7  convolutional  code  v/ith  eight 
level  soft  Viterbi  decoding  was  demonstrated  by  the  simulation. 
The  simulation  results  were  within  0.7  dB  of  theoretical. 

All  computer  simulations  were  done  at  baseband  to  reduce 
simulation  times.  The  Monte  Carlo  error  counting  technique  was 
used  to  estimate  P^.  The  effect  of  increasing  the  samples  per 
bit  in  the  simulation  was  demonstrated  by  the  0.4  dB  improvement 
in  Pb  caused  by  doubling  the  number  of  samples. 
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COMPUTER  SIMULATION  OF  A  GENERAL  PURPOSE  SATELLITE  MODEM 


I.  Introduction 


1. 1  Background 

The  modern  world  of  communications  :,.s  a  rapidly  changing 
environment.  In  the  area  of  satellite  communications,  systems 
using  digital  signal  processing  (DSP)  are  increasingly  found. 

The  CQM-248  modem  is  a  modem  that  uses  DSP  to  implement  many 
functions  [1].  As  communications  system  complexity  has  grown, 
the  analytical  evaluation  of  the  system  performance  has  become 
difficult  if  not  impossible  [2].  To  attempt  to  predict  system 
performance,  engineers  turn  to  computer  simulations  of  a  system 
model.  The  criteria  often  used  for  system  performance 
measurement  is  the  probability  of  bit  error  (Pb). 

1.1.1  Computer  Simulation.  The  availability  of  simulation 
packages  has  made  the  use  of  computers  to  simulate  and  evaluate 
the  performance  of  communication  systems  commonplace  [ 3 ] .  The 
use  of  an  interactive  package  relieves  the  engineer  from  the 
tedious  job  of  computer  programming  and  calculation  of  complex, 
often  repetitive,  calculations. 

The  computer  simulation  is  a  mathematical  model  of  the 
communications  system.  The  computer  performs  calculations  on  the 
independent  inputs  (data  bits)  and  produces  output  (data  bits) 
based  on  the  model.  Each  continuous  time  data  bit  is  represented 
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in  a  scunpled  form  in  the  simulation.  Direct  comparison  of  system 
input  and  output  data  bits  for  differences  can  be  used  to 
estimate  Pi,. 


1.1.2  Monte  Carlo  Technique.  The  Munte  Carlo  technique  is 

an  error  counting  technique  used  to  estimate  [  3  ] .  To  obtain 

an  error  rate  estimate  that  approaches  the  true  Pj,,  many  data 
bits  must  be  compared.  The  rule  of  thumb  is  that  a  minimum  of 
10’'**  bits  must  be  observed  to  estimate  a  true  P^  on  the  order  of 
10"''  [2,3].  The  10’'*^  bits  define  a  confidence  interval  where  the 
estimated  Pt  will  range  from  one  half  to  twice  the  true  P^,.  This 
interval  is  considered  acceptable  [3]. 

1.1.3  COM-248  Modem  Basics .  The  modem  is  designed  to 
implement  phase-shift  keyed  modulation  (PSK)  on  digital  data  [1]. 
The  digital  data  is  scrambled,  differential,  and  convolutional 
encoded  before  wave  shaping  occurs.  The  wave  shaping  is  done  by 
proprietary  application-specific  integrated  circuits  and  is  not 
considered  in  this  thesis.  The  modem  has  both  sequential  and 
Viterbi  decoding  capabilities  for  demodulation  of  the  encoded 
data.  The  modem  is  capable  of  producing  either  binary  PSK  (BPSK) 
or  quadrature  PSK  (QPSK)  forms.  Options  are  also  available  to 
implement  8-PSK  with  trellis  coded  modulation  and  off-set  QPSK 
(OQPSK) . 

1.1.4  Code  Gain .  The  use  of  coding  schemes  in  communication 
systems  has  become  an  accepted  practice  [ 4 ] .  Convolutional 
coding  with  soft  decision  Viterbi  decoding  is  a  standard 
technique  used  over  satellite  communication  channels. 
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Performance  of  Viterbi  decoding  with  eight  level  quantized  soft 
decisions  has  been  shown  to  be  within  0.25  dB  of  theoretical 
predictions  [5].  For  a  rate  1/2  constraint  length  7 
convolutional  code  with  ideal  BPSK  detection,  the  code  gain 
predicted  by  [4]  is  3.8  dB  for  Pb  =  1X10'^. 

1.2  Goal 

The  goal  of  this  thesis  is  to  estimate  the  CQM-248  modem 
performance.  This  will  be  accomplished  through  the  use  of 
computer  modeling  and  simulation.  The  probability  of  bit  error 
in  the  presence  of  additive  white  Gaussian  noise  is  the 
performance  measure.  The  Monte  Carlo  technique  is  used  to 
estimate  the  probability  of  bit  error.  Specifically,  the  Signal 
Processing  Worksystem™  (SPW™)  software  tool  is  used  to  model  the 
communications  system.  To  accomplish  this  goal,  three  separate 
objectives  will  be  pursued.  The  first  objective  is  to  gain 
experience  with  the  SPW™  tool.  The  second  objective  is  to 
develop  and  demonstrate  proper  operation  of  a  Costas  loop 
demodulator.  The  final  objective  is  to  assemble  and  test  a  full 
system  model. 

1.3  Assumptions 

Five  assumptions  affecting  the  modem  simulation  are  made. 

1.  The  effects  of  wave  shaping  on  Pb  will  be  neglected. 

Due  to  the  proprietary  wave  shaping  used  in  the  CQM-248,  no 
attempt  is  made  to  model  wave  shaping. 

2.  Perfect  symbol  synchronization  is  achieved. 
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The  symbol  synchronization  used  by  the  modem  is  a  slope  detection 
technique  that  works  on  shaped  pulses  [ 6 ] ,  and  because  shaped 
pulses  are  not  modeled,  any  attempt  at  modeling  symbol 
synchronization  is  unjustified.  Symbol  synchronization  is 
accomplished  in  the  simulation  by  use  of  simulation  generated 
timing  signals. 

3.  The  simulation  random  data  generator  is  assumed  to  be 
sufficiently  random  to  eliminate  the  scrambler  part  of  the  modem 
from  the  model. 

4.  The  automatic  gain  control  (AGC)  portion  of  the  modem  is 
assumed  to  negligibly  affect  probability  of  bit  error  and  is  not 
included  in  the  model.  It  has  been  shown  that  variation  of  AGC 
gain  of  ±20%  is  tolerated  by  Viterbi  decoding  (soft  decision) 
with  no  significant  performance  degradation  [5]. 

5.  Finally,  no  attempt  will  be  made  to  model  the  modem 
signal  acquisition  processes. 

1.4  Scope 

The  simulation  system  is  used  to  model  two  modems  operating 
over  a  channel  that  disturbs  the  transmitted  signal  with  additive 
white  Gaussian  noise  (AWGN).  All  simulation  is  done  at  baseband 
due  to  the  long  simulation  run  time  involved  at  actual  operating 
frequencies.  Bandpass  systems  can  be  successfully  modeled  with 
their  lowpass  equivalent  model,  thus  reducing  the  simulation 
seunple  frequency  and  computer  run  time  [2] . 

The  modem  configuration  modeled  is  the  differentially 
encoded  BPSK  (DEBPSK)  option.  Data  will  have  both  differential 
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and  convolutional  encoding  applied.  The  convolutional  code  used 
is  the  optimum  rate  1/2  constraint  length  7  code  [4].  The 
demodulation  process  of  the  modem  is  modeled  as  a  second-order 
Costas  loop.  The  probability  of  bit  error  will  be  estimated  for 
four  energy  per  bit  to  noise  power  (Eb/No)  levels.  The  lowest 
Eb/No  level  tested  is  3.5  dB,  based  on  practical  experience  with 
the  modem  [6].  Also,  the  probability  of  bit  error  at  lower 
levels  of  Eb/No  for  the  coded  system  rapidly  approaches  uncoded 
BPSK  performance  [ 4 ] . 

1.5  Approach 

1.5.1  BPSK  Matched  Filter  System.  The  first  objective  is 
to  gain  experience  with  the  SPVP*  tool.  To  do  this,  a  basic  BPSK 
matched  filter  detection  system  will  be  modeled  using  standard 
SPW™  system  blocks.  The  AWGN  added  by  the  channel  will  be  used 
to  generate  bit  errors.  The  Monte  Carlo  method  of  estimating  Pb 
will  be  demonstrated.  The  average  of  ten  simulations,  each  run 
with  a  different  noise  seed,  will  be  found  for  three  Eb/N^  levels. 
The  resulting  data  will  be  plotted  and  compared  to  theoretical  Pb 
for  matched  filter  detection  of  BPSK.  There  are  three  goals  to 
this  portion  of  testing.  First,  the  Monte  Carlo  method  is 
demonstrated.  Second,  operation  of  the  bit  error  rate  counter  is 
verified.  Finally,  the  AWGN  generator  block  is  confirmed  by  its 
ability  to  cause  errors  in  the  system  at  the  expected  rate. 

1.5.2  Costas  Loop.  The  second  objective  is  to  develop  and 
demonstrate  the  proper  performance  of  a  Costas  loop  demodulator. 
The  Costas  loop  is  the  key  part  of  the  demodulation  process.  Two 
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test  systems  will  be  built  to  demonstrate  correct  Costas  loop 
operation.  First,  a  system  will  be  built  to  test  the  loop 
response  to  a  phase  step  input.  The  loop  phase  error  rise  time 
to  peak  phase  deviation  and  magnitude  of  the  deviation  are  the 
performance  criteria.  Theoretical  values  for  rise  time  and  peak 
deviation  of  the  phase  error  are  predicted  by  [ 7 ] .  The  second 
test  of  the  Costas  loop  will  be  its  response  to  a  frequency  step 
input.  This  is  commonly  used  to  represent  a  Doppler  shift  in  the 
carrier  frequency  [ 7 ] .  The  time  to  lock  onto  the  frequency  step 
and  steady-state  phase  error  are  the  criteria  used  for  the  test. 
Again,  [7]  provides  the  theoretical  values. 

1.5.3  Full  System  Testing.  The  final  objective  is  to 
assemble  and  test  a  full  system.  The  convolutional  encoder  block 
configured  to  model  the  modem  will  encode  the  differentially 
encoded  random  data.  The  encoded  data  will  be  BPSK  modulated  at 
baseband  and  AWGN  added  to  simulate  the  channel.  The  Costas  loop 
will  demodulate  the  noisy  data  which  will  then  be  soft  Viterbi 
decoded.  The  error  counter  will  compare  the  decoded  data  to  the 
original  random  data  and  total  the  errors.  A  curve  of  vs  Eb/N^ 
will  be  produced  from  the  measured  data  and  plotted  for 
comparison  to  theory. 

1.6  Overview 

The  thesis  is  organized  in  seven  chapters.  Chapter  II  has 
two  main  parts.  First,  a  hardware  description  of  the  CQM-248 
modem  is  provided.  Then,  an  overview  of  the  software  tool  used 
to  build  the  modem  model  is  included  and  simulation  at  baseband 
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is  discussed.  Chapter  III  covers  theoretical  probability  of  bit 
error  for  BPSK  and  DEBPSK.  Also,  theoretical  code  gain  is  given. 
The  chapter  ends  with  a  section  on  Monte  Carlo  error  estimation. 
Chapter  IV  demonstrates  the  Monte  Carlo  method  applied  to  matched 
filter  detection  of  BPSK.  Two  of  the  goals  in  Chapter  IV  are  to 
test  the  bit  error  counter  and  AWGN  generator  blocks .  Chapter  V 
provides  Costas  loop  theory  and  results  of  testing  the  simulation 
Costas  loop.  Chapter  VI  describes  the  modem  model  testing  and 
results.  Chapter  VII  is  a  summary  of  the  thesis  work  and 
conclusions.  Appendix  A  includes  SPW’'"-generated  system  plots  and 
signal  plots.  Finally,  Appendix  B  contains  data  tables  from  the 
computer  simulation  runs. 
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II .  Eguipmeti't  Descript.ion 


2,1  COM-248  Digital  Satellite  Modem 

The  CQM-248  is  a  digital  PSK  modem  used  in  satellite 
communication.  The  modem  serves  as  a  link  between  the  user's 
baseband  data  terminal  equipment  and  the  intermediate 
frequency  (IF)  of  the  radio  equipment.  It  is  constructed  with  a 
modular  architecture  to  allow  selection  of  operating 
configuration  and  ease  of  maintenance  [1].  A  simplified 
modulator  block  diagram  is  shown  in  Figure  1,  with  the 
demodulator  block  diagram  shown  in  Figure  2 .  The  next  two 
sections  discuss  the  blocks  included  in  the  simulation. 

2.1.1  Modulator  Block  Diagram.  The  modulator  consists  of 
nine  blocks  as  follow:  data  scrambler,  differential  encoder, 
convolutional  encoder,  digital  filter,  digital-to-analog 
converter,  modulator,  IF  synthesizer,  power  control,  and  bandpass 
filter.  The  data  is  first  scrambled  to  reduce  the  chance  of 
transmitting  a  string  of  zeros  or  ones.  For  the  simulation,  a 
random  data  generator  is  used  as  the  input,  so  the  scrambler 
block  is  not  included.  The  modem  uses  differential  encoding  to 
eliminate  the  180  degree  phase  eutibiguity.  Differential  encoding 
is  included  in  the  simulation.  The  next  block  in  the  data  path 
is  the  convolutional  encoder.  The  CQM-248  can  implement  two  code 
rates  with  a  constraint  length  of  either  7  or  9.  The  code  rates 
are  1/2  or  3/4.  Only  the  rate  1/2  constraint  length  7  is 
simulated.  The  generator  functions  of  the  convolutional  code  are 
denoted  as  Go  =  1111001  (binary)  and  Gj  =  1011011  (binary)  [6]. 
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SCRAMBLER 


DIFFERENTIAL 


CONVOLUTIONAL 


Figure  1.  CQM-248  Modulator 


These  functions  are  coinmonly  known  as  the  optimtim  functions  for  a 
rate  1/2  constraint  length  7  code  [8].  A  SPW™  convolution 
encoder  block  was  configured  to  model  this  portion  of  the 
modulator. 

The  next  block  in  the  data  path  is  the  wave  shaping  digital 
filter.  Due  to  the  proprietary  nature  of  the  wave  shaping,  the 
information  needed  to  simulate  this  block  was  not  available. 
Therefore,  it  is  not  simulated.  The  output  from  the  filter  stage 
is  an  in-phase  and  quadrature  signal  that  is  converted  from 
digital  to  analog  and  sent  to  the  QPSK  modulator  block.  The 
digitally  synthesized  carrier  is  modulated  by  the  antipodal  data 
stream.  To  transmit  a  BPSK  signal,  the  same  data  bits  are  sent 
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to  both  the  in-phase  and  quadrature  channel  of  the  modulator. 

For  the  simulation,  a  0  Hz  carrier  is  used  and  only  BPSK 
modulation  is  simulated. 

The  PSK  modulated  signal  power  is  set  by  a  power  control 
block.  Finally,  a  bandpass  filter  bandlimits  the  signal. 

Neither  the  power  control  or  bandpass  filter  are  implemented  in 
the  modem  model. 

2.1.2  Demodulation  Block  Diagram.  Figure  2  is  a  simplified 
demodulator  block  diagram.  For  reference,  a  simplified  block 


filter  recovered 

DATA 


Figure  2.  CQM-248  Demodulator 


diagreun  for  a  Costas  loop  is  shown  in  Figure  3.  Six  modem  blocks 
are  incorporated  into  the  Costas  loop.  The  blocks  included  are 
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Figure  3.  Simplified  Costas  Loop 


the  QPSK  demodulator;  IF  synthesizer;  analog  to  digital 
converters;  digital  filter;  and  demodulator  processor.  Figure  4 
shows  where  the  Costas  loop  fits  into  the  modem  data  flow.  The 
input  filter  is  a  bandpass  filter  capable  of  passing  the  range  of 
intermediate  frequencies  available  from  the  CQM-248.  The  range 
of  frequencies  is  52  MHz  to  88  MHz  or  104  MHz  to  176  MHz  [1]. 

This  filter  is  not  included  in  the  model  because  the  bandwidth  of 
the  digital  filter  accomplishes  the  noise  spectrum  limiting  [6]. 
Also,  the  automatic  gain  control  (AGC)  is  not  included  in  the 
model  because  it  is  shown  by  [4]  to  have  little  effect  on  bit 
error  rate  for  a  ±20%  range  of  AGC. 
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Figure  4.  Demodulator  Block  Diagram  With  Costas  Loop 


The  modem  block  QPSK  demodulation  corresponds  to  the 
in-phase  (I)  arm  and  quadrature  (Q)  arm  mixers  of  the  Costas 
loop.  The  IF  synthesizer  can  be  considered  the  voltage 
controlled  oscillator  (VCO)  of  the  Costas  loop.  The 
analog-to-digital  conversion  of  the  modem  is  not  simulated 
because  the  simulation  is  all  based  on  sampled  data  bits. 

The  digital  filter  in  the  modem  is  a  patented  application- 
specific  integrated  circuit.  The  digital  filter  is  progreuranable 
to  realize  a  variety  of  equalized  filter  shapes  for  data  rates 
from  below  9.6  Kbps  to  above  2.2  Mbps  [1].  The  demodulator 
processor  is  another  application-specific  integrated  circuit. 

The  demodulator  process  performs  both  carrier  phase  detection  and 
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symbol  synchronization,  providing  control  signals  to  the  digital 
synthesizer  and  closing  the  loop  of  the  carrier  tracking  circuit. 
Two  other  control  outputs  come  from  the  demodulator  processor:  a 
digital  control  signal  used  by  the  AGC  circuits,  and  a  symbol 
synchronization  sent  to  the  digital  filter  and  provided  as  a 
modem  output. 

For  the  model,  2-pole  Butterworth  lowpass  filters  were 
selected  for  the  Costas  loop  arm  filters  for  three  reasons.  The 
Butterworth  arm  filter  performance  in  a  Costas  loop  is  well 
documented  in  [9].  The  modem  digital  filter  information  is 
unavailable.  Lastly,  the  SPW™  Costas  loop  block  was  available 
with  Butterworth  arm  filters.  Phase  detection  is  accomplished 
digitally  in  the  modem  and  will  be  modeled  as  a  mixer.  The  loop 
filter  in  the  demodulator  processor  is  a  proportional  plus 
integrator  circuit  that  can  be  approximated  by  a  lead-lag 
filter  [7].  The  digital  synthesizer  is  modeled  as  a  VCO.  A  full 
discussion  of  the  Costas  loop  demodulator  used  in  the  simulation 
is  provided  in  Chapter  V. 

The  data  output  by  the  demodulator  processor  is  a  quantized 
representation  of  the  analog  input.  Both  the  modem  digital 
filter  and  demodulator  processor  affect  the  quantization  level. 
The  final  result  can  be  modeled  as  a  eight  level  quantized  ±1 
signal  with  (-1,  -0.66,  -0.33,  0,  0.33,  0.66,  1)  as  the 
thresholds  [6].  The  quantized  data  bits  are  used  to  perform  soft 
Viterbi  decoding,  and  the  Viterbi  decoded  data  is  the  input  to 
the  differential  decoder  block.  A  SPW™  Viterbi  decoder  block  was 
configured  to  model  the  modem  Viterbi  decoder. 
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A  differential  decoder  was  included  in  the  model.  As 
previously  stated,  a  random  data  generator  was  used  as  an  input 
to  the  system;  therefore,  both  data  scrambler  and  descrambler 
were  not  included  in  the  model. 

2 . 2  Computer  Simulation  Tools 

The  Signal  Processing  WorkSystem™  and  SPW™  are  trademarks  of 
Comdisco  Systems,  Inc.  [10].  The  SPW™  software  package  provides 
an  interactive  computer-aided  tool  for  digital  signal  processing 
simulation.  The  two  major  parts  of  SPW^"  that  were  used  in  this 
thesis  are  the  Block  Diagrcim  Editor  (BDE)  and  the  Signal  Display 
Editor  (SDE).  Two  copies  of  the  software  tool  were  available  for 
use.  One  copy  ran  on  a  Sun- 4™  and  the  other  copy  ran  on  a 
SPARC- 2™.  The  complied  C  code  from  the  simulation  was  also  run 
on  three  SPARC-2™  workstations. 

2.2.1  Block  Diagram  Editor.  The  BDE  is  used  to  create  and 
edit  systems  consisting  of  signal  processing  blocks.  Signal 
processing  blocks  are  stored  in  libraries.  The  graphical 
interface  of  the  BDE  allows  graphic  representations  of  the  blocks 
called  symbols  to  be  connected  to  form  a  system.  Multiple 
symbols  can  be  joined  together  as  a  detail  and  a  new  symbol 
created  that  performs  the  function  of  the  detail.  Signal 
processing  by  a  symbol  can  be  controlled  by  editing  the  symbol 
pareuneters.  Parameters  control  such  operating  functions  as 
sample  frequency  or  bit  rate.  Parameters  can  be  set  to  a  desired 
value  or  calculated  from  other  parameters.  For  example,  a  gain 
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value  in  a  feedback  loop  can  be  set  to  -10  or  calculated  to  be 
-10  from  the  equation:  -(scimple  frequency /bit  rate). 

A  completed  system  with  parameters  set  to  the  desired  values 
or  equations  can  be  interactively  simulated  from  the  BDE.  When 
the  simulation  run  is  started,  the  number  of  samples  to  process 
(iterations)  is  entered  and  the  default  noise  seed  of  1  is  used. 
The  results  of  the  simulation  run  are  stored  in  files  in  two 
ways.  The  signal  sink  block  in  a  system  will  store  numeric 
values  in  a  file  that  can  be  displayed  by  the  SDE.  The  write 
results  block  stores  the  block  instance  number  and  iteration 
number  along  with  the  system  result  to  a  file.  The  write  results 
block  output  can  be  viewed  in  a  text  window  in  the  BDE  or  with  a 
text  editor  on  the  computer. 

2.2.2  Signal  Display  Editor.  The  SDE  software  module  is 
used  to  create,  edit,  display,  and  analyze  signal  waveforms. 
Signals  stored  from  a  simulation  run  can  be  overlaid  to  check 
timing.  Also,  math  functions  can  be  performed  on  the  stored 
signals  to  normalize  results.  The  SDE  allows  all  iterations  to 
be  displayed  at  one  time  or  a  few  iterations  to  be  zoomed  in  on 
to  find  a  peak  value.  Typical  plots  from  the  SDE  are  shown  in 
Appendix  A. 

2 . 3  Simulation  at  Baseband 

In  general,  a  modulated  carrier  can  be  represented  in 
quadrature  form  given  by  [11]  as: 


15 


Sit)  =  ATj  ( t)  cos  (27tfjjt  +  6)  -  ^2  ( t)  sin (27tfpt  +  0) 


(1) 


where 

X,(t)  and  ^^(t)  are  low-pass  processes 
fo  is  the  carrier  frequency 
B  is  the  bandwidth  of  the  low-pass  process 
fo  is  typically  much  greater  than  B,  and 

0  is  the  carrier  phase  offset  also  called  channel  rotation 

The  complex  envelope  form  of  Equation  1  is 

Sit)  =[J:i(t)  +jX2(t)]e^®  (2) 

where  f^  =  0  Hz. 

In  the  simulation,  the  continuous  time  signal  S(t)  is 
processed  as  a  scimpled  signal  S(kAt)  with  a  sampling  frequency 
f,,  *  1/At,  and  k  is  the  particular  sample  number.  In  SPW™,  the 
sample  number  is  referred  to  as  the  iteration  number.  The 
representation  of  S(t)  for  a  given  sampling  frequency  f„  is  given 
as  [ 2  ] : 

SikAt)  =  AikAt)/Pe<  f  (3) 

where 

A(kAt)  represents  amplitude  modulation 
<E>(kAt)  represents  phase  modulation 
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is  the  relative  carrier  frequency 
P  is  the  signal  power 

is  the  sampling  frequency,  and 
k  is  the  iteration  index 

For  BPSK,  S(kAt)  will  have  the  form  [2]: 


5(kAt)  =  A 


(4) 


where 

A(kAt)  is  a  random  ±1  modulating  bit  stream,  and 
0  is  a  constant  channel  rotation 


For  the  baseband  model,  f^  =  0  Hz.  The  sampling  rate  can  now  be 
set  based  on  the  modulating  symbol  rate  and  the  effects  of 
aliasing  on  the  power  spectra  aensity.  As  discussed  in  [2],  the 
f,a  producing  the  best  results  will  be  an  even  integer  8  to  16 
times  the  symbol  rate.  If  less  than  8  scunples  per  symbol  are 
used,  accuracy  is  lost,  and  with  more  than  16  samples  per  symbol, 
simulation  time  may  become  excessive  with  little  gain  in 
accuracy  [2].  For  simulation  of  the  modem,  f„  is  set  to  10 
times  the  bit  rate. 
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2.4  Summary 


This  chapter  provided  an  introduction  to  the  modem  to  be 
modeled.  The  data  flow  through  the  modem  was  discussed,  and  the 
major  SPW*"  blocks  used  in  the  simulation  related  to  modem 
blocks.  A  major  point  was  the  way  the  Costas  loop  incorporates 
blocks  from  the  modem.  Also,  the  fact  that  information  needed  to 
simulate  the  proprietary  wave  shaping  was  not  available  was 
mentioned.  The  Block  Diagram  Editor  and  Signal  Display  Editor 
from  the  software  tool  used  for  the  simulation  were  presented. 
Finally,  computer  simulation  at  baseband  considerations  are 
presented,  and  the  method  used  to  select  the  sampling  frequency 
was  discussed.  A  saunpling  frequency  ten  times  the  bit  rate  will 
be  used  in  simulation. 


18 


This  chapter  presents  the  theoretical  probability  of  bit 


error,  Pf  Theoretical  curves  for  Pt  vs  energy  per  bit  to  noise 
power  (Efa/No)  ratio  are  plotted  for  binary  phase  shift  keying 
(BPSK)  and  differentially  encoded  BPSK  (DEBPSK).  Next,  the 
theoretical  code  gain  for  convolutional  encoding  is  reviewed.  A 
plot  of  Pb  vs  Efc/No  resulting  from  the  convolutional  encoding  is 
provided.  Then,  the  effect  of  imperfect  carrier  reference  on  the 
Pb  curve  is  discussed.  The  chapter  ends  with  the  theory  behind 
the  Monte  Carlo  method  of  Pb  estimation.  The  confidence  interval 
associated  with  the  use  of  the  Monte  Carlo  Pb  estimate  is 
presented. 


Probability  of  Bit  Error 


A  common  measure  of  performance  used  for  comparing  digital 
modulation  methods  is  the  probability  of  symbol  error  (P^)  [8]. 
For  BPSK,  the  probability  of  symbol  error  is  the  same  as  the 
probability  of  bit  error  Pb-  For  coherent  detection  of  BPSK  with 
antipodal  signaling,  the  equation  for  Pb  is  given  in  [8]  as 


where 

Eb  is  the  energy  per  received  bit 
No  is  the  white  noise  power 

Q(’ )  is  the  complementary  error  functions  defined  by 
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(6) 


0(x) 


1 

v'Sii’ 


X 


The  energy  per  bit  is  [8] 


T  r 

=  fsUt)dt  =  f  slit)  dt 

0  0 

For  antipodal  BPSK,  where  S(t)  is  the  bit  cunplitude  ±A,  and 
S?(t)  =  S2(t)  =  S^(t),  Ej,  reduces  to 


T 

Et,  =  ^A^dt  =  A^T 
0 


(8) 


When  suppressed  carrier  recovery  methods  are  used  for  BPSK 
carrier  recovery,  a  ±180*  phase  ambiguity  exists  [8]. 

Differential  encoding  is  a  common  method  used  to  avoid  the  phase 
ambiguity.  The  probability  of  bit  error  for  coherently  detected, 
DBBPSK  is  given  by  [ 8 ] . 


Pb  =  20 


2E^ 


N 


1-0 


2Et, 

N  } 


(9) 


Figure  5  is  a  plot  of  Equation  5  and  Equation  9 .  The  plot 
was  computer  generated  for  1/4  dB  steps  of  the  independent 
variable  Eb/Ng. 
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Eb/N. 


Figure  5.  Theoretical  Probability  of  Bit  Error 


3.2  Theoretical  Code  Gain 

The  effect  of  using  a  convolutional  code  on  the  bit  stream 
before  BPSK  modulation  and  soft  decision  Viterbi  decoding  on  the 
demodulated  bit  stream  is  discussed  in  [4,5].  The  desired  result 
of  coding  is  to  reduce  the  amount  of  signal  energy  required  to 
receive  data  at  a  specified  Pt.  The  amount  of  reduction  of  Eb/N^ 
required  to  achieve  a  given  Pt  is  defined  as  the  code  gain  [4]. 
For  the  rate  1/2,  constraint  length  7  code  with  the  same 
generator  polynomials  used  in  the  modem,  [4,5]  predict  a  coding 
gain  shown  in  Table  1. 
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Table  1.  Pred: 

.cted  Code  Gain 

Pb 

Code  Gain 

1 

o 

3.8  dB 

10-® 

5.5  dB 

lo-’ 

5.8  dB 

Upper  Bound 

7.0  dB 

This  predicted  gain  is  from  simulation  with  equally  spaced 
eight  level  quantization  for  (-1.5,  -1,  -0.5,  0,  0.5,  1.0,  1.5) 
thresholds.  The  path  length  of  the  Viterbi  decoder  was  32. 
According  to  [4,5]  the  eight  level  quantization  suffers  a  loss  of 
less  than  0.25  dB  compared  to  infinitely  fine  quantization. 

Also,  a  path  length  in  the  decoder  of  4  to  5  times  the  code 
constraint  length  is  sufficient  for  negligible  variance  from 
optimum  decoder  performance. 

In  Figure  6,  the  curve  labeled  DEBPSK  with  coding  was  made 
by  subtracting  the  predicted  code  gain  from  the  theoretical  P^, 
curve  generated  from  Equation  9 .  The  theoretical  code  gain  curve 
will  be  compared  to  the  simulation  performance. 

The  theoretical  code  gain  assvunes  perfect  carrier  phase 
synchronization.  In  real  systems  an  imperfect  carrier  phase 
reference  degrades  the  Pt  performance  of  the  system.  A  coded 
system  is  particularly  sensitive  to  tracking  errors  due  to  the 
steepness  of  the  Pt  versus  E^/No  curves  [5].  For  suppressed 
carrier  recovery  by  a  second-order  phase- locked  loop,  the  loop 
signal-to-noise  ratio  (pi)  effect  on  Pj,  was  shown  in  [5]  to 
degrade  the  performance  by  approximately  0.5  dB  for  P^  =  10‘®  and 
Pi  “  13.5  dB.  The  degradation  is  measured  from  the  pi  =  <»  result. 
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For  Pi  =  12  dB  at  Pb  =  10‘®,  the  degradation  was  shown  to  be 
approximately  3  dB,  whereas  for  pi  =  15  dB,  the  degradation  was 
less  than  0.5  dB. 

3.3  Methods  of  Error  Estimation 

A  method  of  error  estimation  is  required  for  many  digital 
systems  because  of  nonlinear  operations  in  the  demodulation 
process  [2].  For  example,  the  quantizing  of  the  sampled  signals 
and  squaring  of  the  signal  for  carrier  recovery  are  both 
nonlinear  process. 

The  demodulation  and  detection  process  reduces  the  waveform 
to  a  number  which  is  compared  to  a  threshold.  The  decision 
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process  can  be  described  in  terms  of  the  probability  density 
functions  (pdf)  fo(V;T)  and  /i(V;x)  of  the  input  voltage  at  the 
Seunpling  instant  t.  The  decision  that  a  "one"  or  a  "zero"  was 
sent  is  made  based  on  the  input  voltage  at  the  sampling  time  x . 
Two  possible  errors  exist.  The  decision  that  a  "one"  was 
received  when  a  "zero"  was  sent  indicating  a  large  positive 
excursion  of  the  received  voltage  from  a  value  of  a  "zero". 
Similarly,  a  sufficient  error  in  the  value  of  a  "one"  will  cause 
the  error  in  decision  resulting  in  the  received  value  being 
declared  a  "zero".  Sample  pdf  plots  are  shown  in  Figure  7. 
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Figure  7.  Example  Probability  Density  Functions  [3] 


In  the  most  general  sense,  fo(V)  does  not  have  to  equal  the 
shape  of  fi(V).  To  simplify  notation,  the  x  dependence  is 
dropped.  The  probability  of  error,  given  that  a  "one"  was  sent, 
is 
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(10) 


Prob [error /one] 


The  probability  of  error,  given  that  a  zero  was  sent,  is 


Prob  [error/zero]  -  Po  =  jfa(v]dv 

Vt 


(11) 


The  average  probability  is  then 

P  -  ^iPi  f^oPo  (12) 

where 

Ki  is  the  a  priori  probability  of  the  symbol  "one" 
is  the  a  priori  probability  of  the  symbol  "zero" 

Five  methods  of  making  an  estimate  of  the  bit  error  rate  are; 

a)  Monte  Carlo  simulation 

b)  modified  Monte  Carlo  simulation 

c )  extreme-value  theory 

d)  tail  extrapolation 

e)  quasi-analytical  (simulation  run  without  noise  combined 
with  analytical  representation  of  noise). 

The  methods  vary  in  their  assumptions  about  and  methods  for 
working  with  the  pdfs.  Because  the  Monte  Carlo  method  makes  no 


25 


a  priori  assumptions  about  fo(V)  and  fi(V)  and  is  the  most  general 
of  the  five  methods,  it  was  selected  for  use  in  the  simulations. 
Also,  the  nonlinearity  introduced  by  the  quantization  in  the  soft 
Viterbi  decoder  makes  analytical  methods  impractical.  For  a  in- 
depth  discussion  of  the  methods,  see  references  [3,11]. 

The  Monte  Carlo  method  relies  on  error  counting  to  estimate 
the  probability  of  bit  error.  Assvime  a  "zero"  was  sent.  The 
probability  of  error  is 


Po  =  ffo(v)dv 

''t 


(13) 


Now  consider  an  error  detector  defined  by 


hjv) 


1,  viVj. 

O,  v<Vj. 


This  allows  Equation  13  to  be  rewritten  as 


(14) 


Po  -  f  *o(v)  fo(v)  dv 


(15) 


Equation  15  can  be  viewed  as  the  expected  value  of  h„. 


=  F[/2„(v)] 


(16) 
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The  probability  of  error  can  be  estimated  with  a  sample  mean  of 

ho# 


(17) 

where 

ho  is  the  error  detector 
The  summation  is  an  error  counter 
Vi=V(ti),  t^;  is  the  instant  of  decision  time 
i  is  the  i  th  symbol. 

In  simple  terms,  N  symbols  are  checked  for  error  and  the 
total  error  t)  is  add  up.  Equation  17  defines  the  Monte  Carlo 
method  [3].  Equation  17  can  be  extended  to  account  for  errors 
when  a  "one"  was  sent,  but  detected  incorrectly.  Letting  t) 
represent  the  total  number  of  errors  out  of  N  bits  observed 
defines  an  unbiased  Monte  Carlo  estimate  of  probability  of  bit 
error  p. 


J1 

N 


(13) 


The  probability  estimate  of  bit  error  p,  approaches  the  true  error 
probability  p  as  the  number  of  bits  observed  N  approaches  oo. 

3.3.1  Confidence  Interval.  The  Monte  Carlo  method  allows 
estimating  the  true  probability  of  error  p  within  a  specified 
confidence  interval.  After  the  confidence  interval  is  selected. 
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the  number  of  bits  that  must  be  processed  through  the  system 
simulation  can  be  found  for  a  desired  true  value  P,.  The 
confidence  level  is  defined  as  [3] 


P[h2  i  P  i  hj  =  (1-0)  ^ 

where 

hi-hj  is  the  confidence  interval 
(1-a)  is  the  confidence  level 

The  true  value  of  the  bit  error  rate  p  will  lie  between  hj 
and  hj  with  probability  (1-a).  Figure  8  shows  the  confidence 
interval  for  three  confidence  levels  (90%,  95%,  99%).  Recalling 
the  rule  of  thumb  from  Chapter  I  that  N  should  be  on  the  order  of 
10/P,,  the  horizontal  axis  of  Figure  8  would  be  entered  at  the 
N  =  lO***'  point.  The  90%  confidence  level  at  the  N  =  10*'**  point 
produces  a  confidence  interval  of  approximately  0.5p  to  2p.  This 
interval  is  considered  acceptable  [3]. 

3.3.2  Estimation  Error.  The  normalized  error  of  the 
estimated  p,  is  defined  as  [12] 

standard  deviation  of 

e - p-^ -  (20) 

The  true  bit  error  probability  in  the  system  is  P,  and  can  be 
calculated  analytically  for  some  cases.  An  estimator  is 
considered  acceptable  if  the  normalized  error  e  is  less  than 
1.0  [12]. 
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If  no  assumptions  are  made  about  the  distribution 
determining  the  E[p„]  or  the  variance  of  p,,  use  of  the  Chebyshev 
inequality  is  required  to  calculate  a  upper  bound  on  the 
probability.  The  Chebyshev  inequality  has  the  form  given  by  [13] 
as 


where 

p,  is  a  random  variable 

is  the  standard  deviation  of  p. 
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H  is  the  mean  =  E[p,] 

t  is  the  number  of  standard  deviation  units  away  from 

For  t  =  10,  the  value  of  a  particular  p„  run  will  be  within  10 
standard  deviation  units  of  E[p„]. 

3«3.3  Effect  Of  Dependent  Errors .  The  assumption  that  error 
events  are  independent  does  not  always  hold  true.  A  common 
example  is  differential  encoding  of  the  bit  stream.  The  decoding 
process  will  tend  to  produce  errors  in  pairs.  According  to  [11], 
the  variance  of  the  estimator  will  be  stretched  by  a  factor  of 
( 1  +  2m)  where  errors  occur  in  bursts  of  ( 1  +  m) .  For 
differential  encoding,  m  =  1  and  the  standard  deviation  of  p,  will 

be  multiplied  by  ^  .  Thus  e  will  be  scaled  by  yj  for 

differential  encoding.  Also,  the  confidence  interval  will  not  be 
as  tight  as  expected  for  the  case  of  independent  errors  in  a 
Monte  Carlo  simulation. 

3.4  Summary 

The  theoretical  curves  of  Pj,  vs  Eb/N^  were  plotted  from 
Equation  5  and  Equation  9  for  ideal  recovery  of  BPSK  and  DEBPSK 
respectively.  Next,  the  code  gain  expected  from  the 
convolutional  code  was  presented  and  plotted  for  the  DEBPSK  case. 
Code  gains  of  from  3.8  dB  to  5.8  dB  are  predicted  for  Pt  ranging 
from  10'^  to  lO"’.  Finally,  the  Monte  Carlo  Pt  estimation 
technique  was  discussed.  The  error  counting  done  in  the 
simulation  will  be  used  to  implement  the  Monte  Carlo  technique. 
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with  the  number  of  bits  observed  on  the  order  of  lO/Pt,  a  90% 
confidence  interval  is  achieved  for  independent  errors. 


31 


IV.  BPSK  Matched  Filter  Test 


As  a  starting  point  for  working  with  SPW™,  a  BPSK  system  was 
assembled  and  tested.  An  important  part  of  this  simulation  is 
verification  of  the  SPW™  noise  generator  block  operation.  Next, 
the  configuration  of  the  bit  error  counter  is  given.  Also, 
timing  considerations  for  the  simulation  are  explained.  Then, 
the  testing  process  used  is  covered.  Finally,  the  simulation 
results  are  presented  and  compared  to  the  theoretical  Pt. 

4.1  BPSK  Modulation  in  SPW^" 

A  simplified  block  diagram  of  the  simulated  BPSK  system  is 
shown  in  Figure  9.  The  SPW™  block  diagram  is  shown  in  Figure  17 
in  Appendix  A.  The  SPW™  BPSK  modulator  block  generates  a 
constellation  that  is  tilted  45®  to  (1,  1)  and  (-1,  -1).  The  use 
of  the  complex  tone  generator  block  and  mixer  block  internal  to 
the  BPSK  modulation  block  forces  energy  onto  both  the  in-phase 
and  quadrature  parts  of  the  complex  signal.  When  the  carrier 
frequency  is  set  to  other  than  baseband,  the  output  of  the  BPSK 
modulation  block  was  found  to  have  ±1  data  changes  on  both 
channels,  or  a  45®  tilt,  instead  of  ±1  data  on  the  in-phase 
channel  and  zero  on  the  quadrature  channel. 

As  developed  in  [14],  the  variance  of  the  noise  oj  depends 
on  the  value  of  T]o  and  the  sampling  rate  f„  of  the  simulation. 

The  simulation  bandwidth  is  shown  in  Figure  10.  The  equation  for 
simulation  bandwidth  is 
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where  f„  is  the  sampling  frequency. 


(22) 


The  variance  oj,  or  power  in  a  zero-mean,  real-valued  Gaussian 
random  process  within  a  simulation  bandwidth  B,  is  calculated  to 
be 


Wo  =  = 


(23) 


where  T|o  is  the  white  noise  power  spectral  density. 
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Figure  10.  Simulation  Bandwidth  of  a  Gaussian  Random 
Process 


The  SPW™  block  white  noise  given  signal- to-noise  ratio  and 
single-sided  bandwidth  produces  complex  Gaussian  white  noise 
samples.  The  noise  variance  is  calculated  from  the  following 
equation  [10]: 

noise-variance  =  (signal_average_power/SNR)*(f„/(4*BW) )  (24) 

where 

E 

SNR  is  defined  as  —  with  E,  the  energy  per  symbol 

■N'o 

BW  is  the  symbol  rate 

The  SPW*"  equation,  with  some  manipulation  and  noting  that  two 
times  the  single  sided  noise_BW  gives  the  2-sided  noise 
bandwidth,  is  the  same  as  Equation  23. 
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4.2  Data  Flow  In  The  System 

The  random  data  block  produces  random  bits,  0  for  logic  low 
and  1  for  logic  high,  at  the  system  bit  rate.  The  BPSK 
modulation  block  produces  a  baseband  complex  signal  because  the 
carrier  parameter  was  set  to  0  Hz.  The  0  to  1  logic  input  is 
changed  to  an  antipodal  ±1  V  bit  stream  internally  before  BPSK 
modulation.  Next,  a  discussion  of  noise  generation  is  required. 

The  SPW™  noise  generation  block  detail  is  shown  in  Figure  18 
in  Appendix  A.  It  takes  the  output  from  the  BPSK  modulation 
block  as  an  input.  The  desired  SNR  is  passed  to  the  noise 
generation  block  as  a  parameter  with  units  of  decibel.  The  white 
noise  given  SNR  and  bandwidth  block  output  starts  as  a  zero  mean 
and  unit  variance  complex  white  noise  block  output.  One  block 
input  and  three  parameters  are  used  in  the  block  calculations  of 
scale  factors  required  to  produce  the  desired  noise  variance  as 
given  in  Equation  24.  The  complex  average  power  of  the  noiseless 
signal  input  to  the  block  is  calculated  and  divided  by  the  SNR 
parauneter.  The  sannpling  frequency  and  single-sided  noise 
bandwidth  pareimeters  are  used  to  generate  the  second  scale 
factor.  The  complex  output  from  the  block  is  the  Gaussian  white 
noise  with  zero  mean  and  variance  adjusted  to  provide  the  desired 
SNR. 

For  the  PSK  matched  filter  system,  the  SNR  specified  for  the 
white  noise  block  is  the  Eb/N^  used  to  calculate  the  theoretical 
Pb.  The  noise  is  then  summed  with  the  complex  BPSK  signal.  This 
step  simulates  AWGN  of  the  channel.  The  noisy  BPSK  signal  is 
sent  to  the  SPW™  block  PSK  matched  filter  demodulator.  A  detail 
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of  the  SPW™  PSK  matched  filter  demodulator  block  is  shown  in 
Figure  19  in  Appendix  A.  This  block  assumes  perfect  symbol  and 
carrier  synchronization.  The  constellation  first  angle  is  input 
as  a  parameter.  The  constellation  first  angle  was  set  to  45®. 

The  PSK  matched  filter  demodulation  block  implements  the  matched 
filters  as  integrate-and-dump  blocks.  The  block  integrates  for 
the  bit  duration  calculated  from  the  parameter  bit  rate.  The 
output  from  the  integrate-and-dump  blocks  is  sent  to  a  PSK 
detector  block.  The  PSK  detector  block  will  compensate  for  the 
amount  of  channel  phase  rotation  specified  by  the  paraimeter 
channel  phase  rotation.  The  unrotated  PSK  signal  is  then  sent  to 
the  PSK  quantize  block  internal  to  the  PSK  detector  block. 

The  PSK  quantize  block.  Figure  20  Appendix  A,  computes  the 
phase  of  the  complex  PSK  signal  at  the  end  of  the  bit  time.  The 
calculated  phase  is  then  compared  to  the  reference  established 
internally  by  the  parameters  PSK  modulation  order  and 
constellation  first  angle.  The  output  of  this  block  is  the 
decoded  complex  bit  stream.  The  real  part  of  this  complex  signal 
goes  to  the  modified  real  error  rate  counter  block. 

The  modifications  to  the  real  error  rate  counter  were  made 
to  generate  the  desired  output  of  number  of  errors  counted  and 
number  of  bits  counted.  The  value  of  initial  error  count  was 
changed  from  0  to  1  because  it  was  found  that  the  error  count  was 
one  low  with  the  first  error  event  generating  a  value  of  0  as 
configured  by  SPW™.  The  number  of  bits  counted  is  available  in 
the  block,  but  not  written  as  an  output.  The  write  block  was 
moved  from  the  probability  of  bit  error  location  to  the  desired 
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niimber  of  bits  counted.  The  modified  block  was  saved,  and  a  new 
symbol  generated  for  use  in  other  systems. 

4.3  Timing  Considerations 

In  order  to  correctly  compare  the  transmitted  bits  to  the 
recovered  bits,  the  transmitted  bits  are  delayed  by  10  iterations 
before  entering  the  bit  error  rate  counter.  The  10  iterations  is 
the  time  it  takes  for  the  integrate  and  dump  block  to  produce  the 
first  valid  data.  The  block  parameter  represents  the  block 
length.  The  error  counter  writes  a  result  at  the  end  of  each 
block.  The  hold  input  on  the  error  counter  block  was  used  to  set 
the  location  in  the  bit  of  the  error  check.  Both  the  transmitted 
and  recovered  bit  streams  are  ±1  V  square  waves  entering  the 
block.  The  midpoint  of  the  bit  was  chosen  as  the  time  to  make 
the  comparison.  The  T_0N  block  produces  a  logic  false  until  the 
simulation  time  exceeds  the  parameter  T_0N.  The  inverted  false 
is  a  logic  true  input  preventing  the  block  from  comparing  data 
before  valid  data  is  available.  After  the  simulation  time 
exceeds  T_0N,  the  T_0N  block  produces  a  logic  true  which  is 
inverted  to  a  logic  low,  allowing  the  position  in  symbol  block  to 
take  over  timing.  The  position  in  symbol  block  produces  a  logic 
true  at  the  specified  time  in  the  symbol.  The  block  parameters 
baud  rate  and  symbol  fraction  are  used  to  generate  the  correct 
timing.  The  position  in  symbol  is  set  to  0.5,  corresponding  to 
the  center  of  the  bit. 
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4.4  Testing  Process 


The  BPSK  matched  filter  detection  system  shown  in  Figure  17 
in  Appendix  A  was  tested  for  three  values  of  SNR  (4,  6,  8  dB). 
For  each  value  of  SNR,  the  theoretical  probability  of  bit  error 
Pb  was  calculated  using  Equation  5.  As  required  by  the  Monte 
Carlo  test  method,  the  number  bits  run  in  the  simulation  was 
10/Pb.  With  ten  samples  per  bit,  the  required  number  of 
iterations  was  100/Pb,  with  an  additional  10  iterations  to  fill 
the  matched  filter.  For  the  simulation,  the  arbitrary  bit  rate 
was  5  bps  and  the  sampling  frequency  was  50  Hz.  This  maintains 
the  desired  10  s2Uiiples  per  bit.  The  ratio  sample  rate  to  bit 
rate  is  the  important  feature,  not  the  individual  values  [2]. 

The  carrier  frequency  was  0  Hz .  The  initial  error  count  was  1 . 
T_0N  was  0.2.  Table  2  summarizes  the  parameters  and  the 
parameter  values. 


Table  2.  Test  Parameters 


Parcimeter 

Name 

Value 

Scimple  frequency 

Sf  req 

50  Hz 

bit  rate  uncoded 

bitrate 

5  bps 

number  of  bits  to  count 

block 

1,000 

signal-to-noise  ratio 

SNR 

3.5  dB 

first  error  event 

1 

time  to  start  error  count 

T_on 

0 . 2  sec 

Table  3  shows  the  three  SNR  values,  block  lengths,  and 
theoretical  probabilities  of  bit  errors. 
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Table  3.  Theoretical 


SNR 

Block  Length 

Number  of  Iterations 

4.0  dB 

1,000 

10,010 

1.25008*10-^ 

6.0  dB 

10,000 

100,010 

2.38829*10-^ 

8.0  dB 

100,000 

1,000,010 

1.90908*10-* 

For  each  simulation,  the  number  of  iterations  was  (10  *  block 
length)  +  10.  The  noise  seeds  were  arbitrarily  chosen.  The  ten 
noise  seeds  are  shown  in  Table  4 . 


Table  4.  Noise  Seeds 


Table  23  through  Table  26  in  Appendix  B  show  the  results  of 
the  testing.  Table  5  shows  the  statistics  of  the  test  results. 
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Table  5.  Statistics  Matched  Filter  Detection  BPSK 


Eb/No 

a 

e 

4.0  dB 

1.34  X  10-* 

2.538  X  10-® 

5.038  X  10'® 

0.403 

6.0  dB 

2.84  X  10*® 

3.804  X  10-’ 

6.168  X  10*® 

0.258 

8.0  dB 

2.43  X  10-* 

3.157  X  10*® 

5.618  X  10-® 

0.294 

The  estimation  error  e  is  less  than  1  for  all  three  levels 
of  SNR.  This  is  an  acceptable  indication  of  the  goodness  of  the 
estimator  [12].  The  mean  value  of  the  ten  simulation  runs  was 
plotted  for  each  E^/No  level.  Additionally,  vertical  lines 
representing  the  value  of  ±o  at  each  E^/No  level  were  plotted. 

The  curve  labeled  measured  data  in  Figure  11  is  within  one 
standard  deviation  of  the  theoretical  curve. 

Although  investigation  of  simulation  time  was  not  a  stated 
goal,  it  did  become  important  at  higher  Eb/N^  levels.  Each  noise 
seed  simulation  took  approximately  10  minutes  at  the  4  dB  Eb/N^ 
level.  As  a  comparison,  over  one  and  a  half  hours  were  required 
to  simulate  10®  bits  with  1.0001X10®  iterations  in  the  interactive 
mode  of  the  Block  Diagram  Editor  (BDE) 

4.6  Summary 

The  simulation  of  BPSK  matched  filter  detection  accomplished 
three  things.  First,  operation  of  the  bit  error  counter  was 
verified.  Second,  the  operation  of  the  white  noise  given  SNR  and 
bandwidth  was  confirmed.  The  verification  of  correct  operation 
is  made  by  comparison  of  the  Pb  vs  Eb/N^  curve  resulting  from  the 
simulation  to  the  theoretical  curve.  The  Pb  curve  generated  from 
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Eb/N. 


Figure  11.  Probability  of  Bit  Error 


the  measured  data  was  within  one  standard  deviation  of 
theoretical.  Finally,  the  estimator  error  e  was  shown  to  be  less 
than  one. 
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V.  Costas  Loop 


This  chapter  begins  with  a  discussion  of  the  use  of  the 
Costas  Loop  as  a  demodulator.  Next,  the  mathematical  equations 
defining  Costas  loop  operation  are  covered.  Then,  the  equations 
necessary  for  calculation  of  simulation  parameters  are  presented. 
The  theoretical  portion  of  this  chapter  ends  with  a  table  of 
parauneter  values  initially  used  for  the  simulations. 

In  the  second  half  of  this  chapter,  the  SPW™  Costas  loop 
block  operation  is  covered.  The  simulation  test  systems  used  to 
find  the  Costas  loop  response  to  a  phase  or  frequency  step  are 
discussed.  Next,  the  Costas  loop  block  modifications  are 

discussed.  Finally,  the  modified  Costas  loop  simulation  testing 
is  presented  and  the  results  compared  to  theory. 

5.1  Costas  Loop  Demodulation  Theory 

The  Costas  loop  or,  I-Q  loop,  is  a  method  of  suppressed 
carrier  recovery.  After  initial  carrier  acquisition,  the  Costas 
loop  provides  a  demodulated  bit  stream  on  the  I  arm,  as  shown  in 
Figure  12.  The  noise  spectrum  shaping  of  the  Costas  loop  arm 
filters  can  be  designed  to  be  the  same  as  the  low-pass  equivalent 
of  the  input  filter  Hi(S)  shown  in  Figure  13.  When  this 
condition  is  met,  the  performance  of  the  Costas  loop  is  identical 
to  the  squaring  loop  [15].  The  squaring  loop  removes  the 
modulation  by  squaring  the  received  signal.  The  phase-locked 
loop  (PLL)  then  tracks  the  double-frequency  term  from  the 
squarer.  After  being  divided  by  two,  the  recovered  carrier  is 
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Figure  12.  Reference  Costs  Loop 


available  as  a  coherent  reference  for  demodulation. 

As  developed  in  [15],  the  Costas  loop  with  multiplier 
gain  ^3^  =  i  and  with  signals  labeled  as  in  Figure  12  has  an  input 
given  by 


X(t)  =  ■/2Sm{C)  sin^  {t)  +  n{t) 


where 

4>(t)  «  +  e(t) 

m{t)  is  the  signal  modulation 
n(t)  is  the  additive  noise 


(25) 
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This  received  signal  is  mixed  with  the  output  of  the  voltage 
controlled  oscillator  (VCO).  The  resulting  in-phase  signal  is 
given  as 


Z^it)  =  X(t)[2(y3Qcos*(t)] 

where 

4f(t)  =  a)„t  Q(t)  is  the  estimate  of  ^(t) 
k(t)  is  the  VCO  output  power 

Substituting  Equation  25  into  Equation  26  and  filtering  the 
double  frequency  terms  yields: 
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y^(t)  =  -  w;,(t)]sin4>(t) 

+  ^2k^Nf.(t)  cos<|)  ( t) 


(27) 


Similarly,  for  Jf(t)  mixed  with  the  quadrature  part  of  the  VCO 
output,  Z,{t)  yields: 

Yg{t)  =  y/5^[v/Sm{ t)  -  Ng ( t) ]sin4) ( t) 

-  y]2k^N^ ( t)  cos(|)  ( t)  ( 2 ® ) 

where 

(t.(t)  -0-0  is  the  phase  error  estimate 

The  phase  detector  mixer  multiplies  the  Y,(t)  and  Y^(t)  outputs 
from  the  arm  filters,  producing  the  phase  error  signal  £(t). 

e(t)  =  kj^[[/sm(t)  -  Ng(t)f  -  Nl{t)}sin2^(t:) 

+  2k^N^(.t)ysw{t)  -  t)]cos24>(  t) 


Finally,  the  loop  filter  closes  the  tracking  loop  by  providing  a 
filtered  version  of  the  e(t)  as  the  control  voltage  input  to  the 
VCO. 

Four  things  are  noteworthy  in  the  previous  equations.  For 
digital  modulation,  m(t)  is  a  ±1  pulse  train;  therefore, 
m*(t)  =  1.  The  loop  is  operating  on  the  square  of  the  noise. 
Also,  twice  the  phase  error  is  being  tracked.  Finally  the 
magnitude  of  e(t)  depends  on  the  power  of  the  received  signal  S. 
The  automatic  gain  control  in  an  actual  receiver  keeps  S 
constant.  The  24»  and  squared  noise  term  introduce  a  squaring 
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loss  term  into  the  loop  bandwidth  equation.  The  loop  noise 
bandwidth  and  loop  filter  transfer  function  F(s)  determine  the 
performance  of  the  Costas  loop.  Because  the  heart  of  the  Costas 
loop  is  a  PLL,  equations  for  a  PLL  were  used  to  determine 
parameters  such  as  VCO  gain  and  loop  filter  natural  frequency  in 
the  model. 

Reference  [16]  documents  a  practical  method  of  designing  a 
PLL  once  the  desired  lock  range  is  chosen.  For  the  model,  the 
known  tracking  parameter  was  the  pull-out  range  cOp^.  The  pull-out 
range  is  the  dynamic  limit  of  stable  operation  once  lock  has  been 
established.  As  defined  in  [16] 


»  1.8q„{C  +  1)  poj 

where 

(ot,  is  the  natural  frequency  of  the  PLL 
X,  is  the  damping  factor 

The  lock-in  range  is  defined  as  [16] 

Ao),  =  2Cio„  (31) 

Given  AWp,,  and  C/  the  loop  natural  frequency  can  be  calculated 
using  Equation  30.  Once  (o„  is  known,  a  value  for  AcDl  can  be 
calculated  using  Equation  31. 

Reference  [16]  assumes  the  values  of  K,,  and  K<j  will  be  found 
in  a  specification  for  the  PLL  device.  As  was  done  by  [16],  for 
simulation  on  a  computer: 


46 


KoKd  =  10  w. 


(32) 


where 

Ko  is  the  VCO  gain  constant 
Kj  is  the  phase  detector  gain 


This  relation  was  initially  used  to  meet  the  requirement  that 
KoKj  »  Assuming  Kj  in  the  simulation  is  1,  a  value  for  can 

be  calculated  from  Equation  32.  The  modem  design  expects 
Afpo  =  2JtAa)po  to  be  20  Hz  or  less  [6].  Equation  30  and  Equation  31 
were  used  to  calculate  theoretical  values  of  f„  =  6.5087382  Hz 
and  Aft  =  9.204  Hz  respectively.  Equation  32  yields  a 
theoretical  of  65.087382  Hz/V.  The  loop  filter 
characteristics,  VCO  gain  K^,  and  phase  detector  gain  determine 
(Do.  To  model  the  modem  loop  filter,  a  passive  lead-lag  filter  is 
used.  For  and  Tj  »  l/K^Kd,  the  passive  filter  approaches 

the  perfect  integrator  form  of  a  proportional  plus  integrator 
loop  filter  [ 7 ] . 

A  passive  lead-lag  loop  filter  has  the  transfer  function 
given  by  [7]  as 


F{s) 


SX2  +  1 
STj  +  1 


where  Xj  and  Tj  are  the  time  constants  of  the  filter. 


(33) 


For  a  passive  loop  filter  of  the  form  in  Equation  33,  [16]  gives 
the  equations  for  and  ^  as 
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(34) 


C  = 


Tl  +  Tj 


(35) 


The  term  K^Kd  is  called  the  loop  gain.  According  to  [16],  the 
loop  is  a  high  gain  loop  if  If  a  PLL  with  a  passive 

loop  filter  is  a  high-gain  loop,  the  phase  transfer  function  can 
be  approximated  as  [7,16] 


H(s) 


+  2i^Oj,S  +  On 


(36) 


This  is  the  same  transfer  function  obtained  if  the  PLL  loop 
filter  is  an  ideal  proportional  plus  integrator  filter  [16]. 

As  was  done  by  [17],  the  Costas  loop  model  operation  was 
verified  using  a  noiseless  simulation.  The  transient  phase  error 
to  step  changes  in  phase  and  frequency  were  the  test  parameters. 
For  a  step  in  phase,  [7,17]  predicted  a  phase  error  response 
where  time  to  peak  deviation  was  used  as  the  criteria  for  proper 
performance.  For  a  step  in  frequency,  the  time  for  the  loop  to 
lock  and  steady-state  phase  error  were  used  as  the  criteria  for 
correct  performance.  The  phase  step  must  be  small  enough  for  the 
loop  to  remain  in  the  linear  region  of  operation  [ 7 ] .  The  small 
angle  approximation  sin0  »  0  for  0  s  30®  is  commonly  used  to 
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determine  the  linear  range.  The  time  to  the  peak  phase  error  for 
the  Costas  loop  using  a  second-order  loop  with  ^  =  0.707,  a  phase 
step  of  20®,  and  f^  =  6.5087  Hz  is  0.0543  sec.  For  testing  the 
loop  response  to  a  frequency  step,  the  frequency  step  must  be 
less  than  or  equal  to  AtOj,.  If  the  frequency  step  exceeds  Acoi,,  the 
loop  will  unlock  and  a  pull-in  process  will  occur  instead  of  the 
desired  lock-in  process.  The  time  to  pull  in  of  a  PLL  is 
typically  1000  times  larger  than  the  time  to  lock  [16].  The  time 
for  a  PLL  to  lock  (Tj,)  onto  a  step  in  frequency,  according  to 
[ 7 ] ,  can  be  approximated  by 


1 


(37) 


For  f„  =  6.5087  Hz,  Tj,  "  0.1536  sec.  The  steady-state  phase  error 
for  a  second-order  loop  using  a  passive  loop  filter  is  given  as 


e8(®®) 


KoKa 


where  Ao)  is  the  frequency  step  in  radians. 


(38) 


For  the  simulation  with  Aco  =  Acoj,  =  2n:Afi,  =  9.2047  Hz  and  = 
65.087  Hz/V,  0e('»)  =  0.1414  radian.  Equation  34  for  and 
Equation  35  for  ^  were  solved  for  expressions  for  and  Tj.  The 
resulting  expressions  are 
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X  =  2C  _ 

"  «„  (40) 

The  noise  bandwidth  and  signal- to-noise  ratio  of  the  loop 
were  the  last  considerations  of  the  theoretical  operation  of  the 
loop.  In  a  Costas  loop  with  two-pole  Butterworth  arm  filters, 
the  squaring  loss  is  approximately  -3  dB  [9].  From  [9]  the 
equivalent  loop  signal-to-noise  ratio  is 


(41) 


with 


(42) 


where 

SNRi  is  the  signal-to-noise  ratio  into  the  loop 
Bi  is  the  arm  filter  bandwidth 

Bi,  is  the  single-sided  loop  noise  bandwidth  of  the  linear 
loop,  given  by: 
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(44) 


Examination  of  Equation  41  for  pi  with  substitution  of 
variaJales  by  Equation  42  and  Equation  43  reveals  that  it  changes 
by  an  amount  directly  proportional  to  SNRi.  Therefore,  a 
constant  can  not  be  maintained  for  the  Costas  loop  suppressed 
carrier  recovery  loop  [ 9 ] .  For  the  theoretical  values  of  Tj  and 
^  given  in  Table  6,  Bj,  =  14.9887  Hz.  The  arm  filter  bandwidth  = 
9600  Hz,  and  SNR^  ranges  from  3.5  dB  to  5.5  dB.  The  theoretical 
range  for  pi  is  approximately  19.5  dB  to  21.5  dB. 

5.2  Summary  of  Theoretical  Operation  Costas  Loop 

Working  with  the  assumption  that  the  modem  acquisition 
circuits  will  perform  the  pull-in  operation,  the  pull-out  range 
for  a  locked  loop  was  the  starting  point  for  the  calculation  of 
loop  parameters.  The  parameters  of  interest  for  the  simulation 
were  loop  damping  loop  natural  frequency  f„,  and  VCO  gain  K^. 
Given  Afpo  =  20  Hz  and  C  “  0.707,  Table  6  shows  the  theoretical 
values  for  loop  operation. 

5 . 3  The  SPW™  Costas  Loop  Block 

The  Costas  loop  block  as  implemented  by  the  SPW™  software 
requires  the  user  to  specify  seven  parameters.  Figure  21  in 
Appendix  A  is  the  unmodified  Costas  loop  block.  Scimpling 
frequency,  natural  frequency  in  Hz,  damping  factor,  and  VCO 
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Table  6.  Theoretical  Costas  Loop  Parameter  Values 


Parameter 

Value 

C 

/2/2 

20  Hz 

fn 

6.5087382  Hz 

Tl 

0.1536396  sec 

Af. 

9.2047  Hz 

KoK, 

65.087382  Hz/V 

0.2123892  sec 

0.0321358  sec 

constant  in  Hz/V  are  all  used  to  control  the  operation  of  the 
loop  filter.  The  loop  filter  is  a  passive  lead- lag  filter.  The 
VCO  quiescent  frequency  and  the  VCO  gain  are  passed  to  VCO  block 
and  used  to  control  its  operation.  The  arm  filter  order  and  arm 
filter  bandwidth  in  Hz  are  passed  to  the  complex  Butterworth 
lowpass  infinite  impulse  response  filter  block.  A  complex 
multiplier  block  is  used  to  mix  the  complex  input.  The  gain  of 
this  block  is  unity.  After  the  complex  output  from  the 
Butterworth  filter  is  separated  into  real  and  imaginary  parts,  a 
real  unity  gain  multiplier  performs  the  phase  detection  process. 
The  complex  conjugate  bock  produces  the  desired  e’^®  signal  from 
the  VCO  output.  The  multiplied  by  e^®  input  removes  the 
channel  rotation  shown  in  Equation  2.  Finally,  the  one  iteration 
hold  block  between  the  lead-lag  filter  and  the  VCO  allows  time 
for  calculations  and  outputs  to  be  generated  in  the  feedback 
path. 
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The  first  modification  made  to  the  Costas  loop  block  was  the 
addition  of  output  ports  for  the  in-phase  and  quadrature  parts  of 
the  Butterworth  filtered  signal.  These  outputs  correspond  to  the 
real-arm  and  imaginary-arm  signals,  as  shown  in  Figure  22  in 
Appendix  A.  Also  added  was  an  output  port  allowing  access  to  the 
phase  error  from  the  phase  detector.  The  block  with  the  above 
modifications  was  saved  and  a  new  symbol  created  to  represent  the 
new  block  as  shown  in  Figure  23  in  Appendix  A. 

5.3.1  Costas  Loop  Testing.  For  testing  the  Costas  loop,  two 
systems  were  built.  Figure  24  in  Appendix  A  shows  the  phase  step 
test  system  with  final  parameter  values.  The  phase  step  system 
multiplies  the  BPSK  signal  by  a  constant  e^®  term.  The  frequency 
step  test  system  simulates  a  frequency  step  by  offsetting  the  VCO 
center  frequency  by  an  amount  in  Hertz  equal  to  the  desired 
frequency  step.  Figure  25  in  Appendix  A  shows  the  frequency  step 
test  system.  The  operation  and  test  results  of  these  two  systems 
are  given  in  the  following  sections. 

5.3.2  Phase  Step  Test  System  and  Results .  The  BPSK  signal 
to  be  multiplied  by  the  e^®  term  was  generated  with  the  BPSK 
modulation  block  fed  by  a  random  data  block.  The  carrier 
frequency  in  the  BPSK  block  was  set  to  0  Hz,  and  the  bit  rate  of 
the  random  data  block  was  set  to  2400  bps.  Using  the  rule  of  10 
samples  per  bit  in  the  simulation,  the  sample  frequency  was  set 
to  24000  Hz.  The  constant  block  value  was  also  set  to  24000  and 
connected  to  the  signal  sink  blocks  sampling  frequency  input 
(s_freq).  The  arm  filter  bandwidth  of  the  two-pole  Butterworth 
filter  in  the  Costas  loop  was  set  to  4800  Hz.  Initial  testing 
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was  done  with  the  theoretical  values  for  loop  frequency  and  VCO 
gain.  The  parameters  were  set  to  6.5087  Hz  and  65.087  Hz/v. 
Table  7  summarizes  parameter  settings  for  the  test. 

Table  7.  Phase  or  Frequency  Step  Test  Parameters 


Pareuneter 

Name 

Value 

sample  frequency 

sf  req 

24000  Hz 

bit  rate 

bitrate 

2400  bps 

loop  frequency 

loopfreq 

6.5087382.0  Hz 

VCO  gain 

VCO_gain 

65.087382.0  Hz/V 

arm  filter  bandwidth 

4800.0  Hz 

step  for  phase  test 

value 

to 

o 

• 

o 

o 

step  for  frequency  test 

VCO_FREQ 

9.2047  Hz 

The  phase  step  0  in  radians  was  produced  by  a  constant  block  with 
value  set  equal  to  the  phase  step  in  degrees  fed  through  a 
conversion  block  degree  to  radian.  The  radian  value  of  the  phase 
step  was  fed  as  the  x  input  to  an  e^*  block  thus  producing  the 
desired  e^®  term.  The  VCO  output  of  the  Costas  loop  is  taken 
before  the  complex  conjugate  block  as  shown  in  Figure  22  in 
Appendix  A.  The  complex  phase  block  converts  the  VCO  output  to  a 
phase.  Correct  Costas  loop  operation  will  produce  a  signal 
stored  in  the  signal  sink  labeled  billsigs/phase_out  that  after  a 
transient  response  matches  the  phase  step.  Table  8  lists  other 
signal  sink  block  names  and  what  they  recorded. 

Correct  operation  of  the  Costas  loop  would  be  verified  by 
the  phase  error  transient  and  steady-state  response.  The  stored 
signals  phase_in  and  phase_out  were  used  to  calculate  the 
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Table  8.  Phase  Step  Test  Signal  Sinks 


Sink  Name 

Contents 

bpsk_out 

modulated  signal 

phase_in 

phase  step  signal  in  radians 

I_ARM 

in-phase  arm  of  Costas  loop  after 
Butterworth  filter 

Q_ARM 

quadrature  arm  of  Costas  loop 
after  Butterworth  filter 

phase_error 

phase  error  from  phase  detector 

VCO_control 

Lead-Lag  filter  output 

VCO_out 

voltage  controlled  oscillator 
output 

phase_out 

phase  of  voltage  controlled 
oscillator  output 

normalized  phase  error.  The  math  functions  in  the  Single  Display 
Editor  part  of  SPW™  were  used  to  perform  the  required 
calculations  on  the  signals. 

The  first  test  runs  of  the  test  system  failed  to  produce  the 
expected  rise  time  to  peak  phase  error.  The  phase  step  size  was 
varied  over  a  range  of  5®  to  20®  with  little  change  in  the  loop 
response.  The  phase  of  the  VCO  output  would  slowly  rise  to  a 
peak  and  oscillate  until  eventually  settling  to  the  phase  step 
value.  This  response  was  similar  to  loops  with  values  of  0.3 
or  less  as  shown  by  [7,16]. 

5.3.3  Frequency  Step  Test  System  and  Results .  The  system 
used  for  testing  the  Costas  loop  response  to  a  step  in  frequency 
is  shown  in  Figure  25  in  Appendix  A.  The  BPSK  modulated  signal 
was  generated  as  a  baseband  signal.  The  frequency  step  was 
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modeled  as  a  frequency  off-set  of  the  VCO  quiescent  frequency. 

The  parcuneter  VCO_FREQ  was  used  to  set  the  VCO  frequency  to  the 
required  step  size.  The  complex  phase  block  was  used  to 
calculate  the  phase  of  the  Costas  loop  VCO  output  during  the 
lock-in  response.  Other  system  blocks  operate  as  described  in 
Section  5.3.2  on  the  phase  step  test  system.  Table  7  summarizes 
pareuneter  settings  for  the  test. 

The  complex  tone  block  was  set  to  generate  a  tone  at  the 
theoretical  loop  natural  frequency  f„.  Comparison  of  this  signal 
to  the  VCO  output  signal  demonstrates  the  lock-in  within  one 
cycle  of  f„.  Table  9  lists  the  signal  sink  names  and  what  was 
stored  in  them. 

Correct  operation  of  the  Costas  loop  would  be  verified  by 
the  time  to  lock  and  steady-state  phase  error  value.  The  test 
system  runs  failed  to  produce  the  expected  results.  For  a  range 
of  frequency  steps  of  3  Hz  to  9.2047  Hz,  the  Costas  loop  appeared 
to  perform  a  pull-in  process  over  time  periods  much  larger  than 
Tf  The  description  given  in  [16]  of  a  pull-in  process  matched 
the  Costas  loop  response.  However,  with  the  Ao)  less  than  Au)l,  a 
lock-in  process  should  have  been  observed. 

5 . 4  Debugging  the  Costas  Loop. 

The  theoretical  calculations  were  checked,  and  the  test 
system  construction  was  thoroughly  checked.  With  the  knowledge 
that  the  Costas  loop  internal  parameter  Xj  and  Xj  were  being 
calculated  by  the  software,  the  equations  used  for  the 
calculations  were  found  and  were 
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Table  9 .  Frequency  Step  Test  Signal  Sinks 


Sink  Name 

Contents 

bpsk_out 

modulated  signal 

loop_f req 

tone  generator  output  at  f„ 

I_ARM 

in-phase  arm  of  Costas  loop  after 
Butterworth  filter 

Q_ARM 

quadrature  arm  of  Costas  loop 
after  Butterworth  filter 

phase_error 

phase  error  from  phase  detector 

VCO_control 

Lead-Lag  filter  output 

VCO_out 

voltage  controlled  oscillator 
output 

phase_out 

phase  of  voltage  controlled 
oscillator  output 

_  VCO  constant 

(|)d  (45) 


,  _  Damping  factor 

•  *3  — 


2 71  VCO  constant 


where 


VCO_constant  is  the  VCO  gain  in  Hz/V 
dainping_f actor  is  the  loop  filter  damping 
f„  is  the  loop  natural  frequency  in  Hz 


(46) 


Comparison  of  these  equations  to  the  theoretical  Equation  39  and 
Equation  40  as  presently  written  is  difficult.  Rewriting 
Equation  39  and  Equation  40  to  change  the  units  from  radian  to 
hertz  produces 
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(47) 


t  -  -±-  .  1 

"  2nfl  2iiKj:a 

t  =  ^  -  1 

=*  2Tif„  2itK^Ka  (48) 

where 

fa  is  the  loop  natural  frequency  in  Hz 

KoKj  is  the  product  of  the  VCO  gain  and  detector  gain  in 
Hz/V 

^  is  the  loop  filter  damping 

With  the  assumption  that  the  detector  gain  is  unity, 

Equations  46  and  48  for  are  identical.  Equations  45  and  47 
for  Xj  are  noticeably  different.  Even  with  the  assumption  that 
Xi»X2,  allowing  the  last  two  terms  of  Equation  47  to  be 
neglected,  the  remaining  equations  do  not  agree.  As  a  test,  the 
theoretical  value  was  written  in  as  the  parameter  value  for  Xj  in 
the  lead-lag  filter  block.  This  effectively  deleted  the  software 
equation  that  had  been  calculating  Xj.  The  Costas  loop  response 
to  the  phase  step  with  the  theoretical  value  of  Xj  in  place  of 
the  software  calculated  value  agreed  with  the  theoretical  rise 
time. 

To  correct  the  Costas  loop  equation  used  to  calculate  Xj, 
the  expression  for  Xj  in  the  lead-lag  filter  was  modified  to 
Implement  Equation  47.  The  modified  expression  for  the  Xi 
pareuneter  was 
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271 


^  _  vco 

■Vi  - 


tant  _  damping  factor  ^ 


T^fr, 


2  IT  VCO_cons  tan  t  (49) 


where 

VCO_constant  is  the  VCO  gain  in  Hz/V 
f„  is  the  loop  natural  frequency  in  Hz 

The  modified  Costas  loop  was  saved  as  a  new  block.  The  phase  and 
frequency  step  test  systems  were  updated  to  use  the  modified 
Costas  loop  and  testing  of  the  loop  completed. 

5.5  Summary  of  Costas  Loop  Testing 

5.5.1  Phase  Step.  For  a  20®  phase  step  with  the  theoretical 
values  for  the  loop  parameters,  Table  10  shows  the  test  results. 
Representative  outputs  of  correct  operation  are  shown  in 
Appendix  A,  Figures  26  -  29. 


Table  10.  Costas  Loop  Test  Result  For  VCO  Gain  65.087382  Hz/V 
and  a  Phase  Step  Input 


Parameter 

Theoretical 

Measured 

rise  time 

0.0543  sec 

0.0544  sec 

normalized  peak  deviation 

-0.208 

-0.167 

steady_state  error 

O 

O 

» 

o 

o 

« 

o 

o 

The  assvimption  of  a  high  gain  loop  in  [  7  ]  did  not  give  a 
value  for  how  much  greater  K^Kd  must  be  than  a)„.  For  the  data 
presented  by  [16],  KoKj  =  10u)„  was  used.  This  also  was  used  for 
the  theoretical  calculations  used  to  calculate  the  VCO  gain 
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pareuneter  resulting  in  Table  10  results.  As  an  additional  test, 
the  simulation  was  repeated  for  the  same  parameters  except  with 
double  the  VCO  gain.  Table  11  shows  the  results  of  using 
130.17476  Hz/V  for  the  VCO  gain. 

Table  11.  Costas  Loop  Test  Result  For  VCO  Gain  130.17476  Hz/V 
and  a  Phase  Step  Input 


Parameter 

Theoretical 

Measured 

rise  time 

0.0543  sec 

0.0559  sec 

normalized  peak  deviation 

-0.208 

-0.192 

steady_state  error 

0 

O 

• 

o 

o 

• 

o 

o 

Rise  time  changed  from  within  0.18%  to  within  2.9%  of 
theoretical.  Peak  deviation  changed  from  within  19.7%  to  within 
7.7%  of  theoretical.  Based  on  the  above  test  results,  the  VCO 
gain  used  for  future  simulations  was  set  to  130.17476  Hz/V. 

5.5.2  Frequency  Step.  The  results  of  the  testing  the 
frequency  step  test  system  are  shown  in  Table  12 .  Testing  was 
done  with  parameter  values  given  in  Table  7  with  the  exception  of 
the  VCO  gain  which  was  set  to  130.17476  Hz/V.  Typical  outputs  of 
correct  operation  are  shown  in  Figures  30  -  33  in  Appendix  A. 

Table  12 .  Costas  Loop  Frequency  Step  Test  Result 


Pareuneter 

Theoretical 

Measured 

time  to  lock 

0.1536  sec 

0.127  sec 

8teady_state  phase  error 

0.0707  radian 

0.0745  radian 
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The  steady-state  phase  error  is  within  5.5%  of  theoretical. 
Figures  30  and  31  in  Appendix  A  show  the  lock-in  process. 

Figure  30  in  Appendix  A  is  the  complex  VCO  output  and  a  tone 
generator  output  at  the  loop  natural  frequency.  The  VCO  center 
frequency  was  offset  by  9.2047  Hz  from  the  0  Hz  carrier.  The  VCO 
output  clearly  shifts  to  a  +1  on  the  real  and  0  on  imaginary 
trace  within  one  cycle  of  the  signal  representing  the  loop 
natural  frequency.  The  VCO  output  then  settles  to  the  steady- 
state  error  value. 

5 . 6  Summary 

This  chapter  covered  the  theoretical  development  of  the 
Costas  loop  parameters.  The  SPW™  block  was  tested  with  the 
theoretical  parameter  values,  and  retested  with  modifications  to 
the  parameter  equations  used  to  calculate  Tj  and  t^.  The  Costas 
loop  with  the  modified  equation  for  Tj  produced  a  simulated 
response  to  a  phase  step  that  was  within  2.9%  of  the  theoretical 
rise  time  and  within  7.7%  of  the  theoretical  peak  deviation.  The 
modified  loop  produced  a  simulated  response  to  a  frequency  step 
that  was  within  5.5%  of  the  theoretical  steady-state  phase  error 
and  lock-in  was  within  one  cycle  of  the  loop  natural  frequency, 
as  predicted  by  theory.  The  VCO  gain  parameter  value  was  doubled 
from  65.087382.0  Hz/V  to  130.17476  Hz/V.  A  trade  was  made  of 
slower  rise  time  for  a  better  peak  deviation  for  the  phase  step 
response.  Rise  time  changed  from  within  0.18%  to  within  2.9%  of 
theory,  whereas  peak  deviation  changed  from  within  19.7%  to 
within  7.7%  of  theory.  Without  the  discovery  and  correction  of 
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the  problem  with  the  SPW™  Costas  loop  equation  for  the  Xj 
parauneter,  the  simulation  of  the  modem  would  not  have  worked 
properly. 
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VI .  Coded  BPSK  System  Simulation 


This  chapter  explains  the  configuration  of  the  SPW™  system 
used  to  model  the  coded  BPSK  system.  First,  the  configuration 
and  testing  of  the  blocks  used  to  encode  the  random  data  are 
explained.  Then,  a  description  of  the  completed  system  used  in 
the  simulation  is  given.  Finally,  the  results  of  the  simulation 
are  presented. 

6. 1  Encoder  and  Decoder  Block  Descriptions 

With  the  Costas  loop  modification  and  testing  completed,  two 
other  systems  remained  to  be  implemented  in  SPW™.  The 
differential  encoder  and  decoder  blocks  were  built  and  tested. 
Also,  the  convolutional  encoder  and  Viterbi  decoder  blocks  were 
configured  to  model  the  modem.  Configuration  and  testing  of  the 
encoder  and  decoder  blocks  are  the  topics  of  the  next  sections. 

6.1.1  Differential  Encoder  and  Decoder.  The  differential 
encoder  and  decoder  blocks  were  built  from  standard  SPW™  blocks. 
Figure  34  and  Figure  35  in  Appendix  A  show  the  encoder  system  and 
symbol  respectively.  The  encoder  takes  in  0  for  a  logic  low  and 
a  1  for  a  logic  high.  The  delay  block  parameter  is  set  to  the 
number  of  samples  in  a  bit.  The  block  was  tested  by  an  input  of 
a  know  bit  stream  producing  a  correctly  coded  output.  The 
differential  decoder  system  and  symbol  are  shown  in  Figure  36  and 
Figure  37  in  Appendix  A.  The  input  to  this  block  must  also  be  0 
for  a  logic  low  and  1  for  a  logic  high.  The  delay  parameter  must 
be  set  to  the  number  of  samples  per  bit.  The  block  operation  was 
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verified  by  its  correct  decoding  of  the  differential  encoder 
block  output. 

6.1.2  Convolutional  Encoder  Configuration .  The  SPW™ 
convolutional  encoder  block  is  well-documented  in  [10].  The 
block  can  be  programmed  to  generate  many  different  convolutional 
codes.  The  block  has  the  following  limits  [10]: 

K  the  number  of  data  bits  input  at  a  time  must  be  <  10 
N  the  number  of  coded  bits  for  each  K  bits  must  be  <  10 
L  the  code  constraint  length  must  be  <  10 

also 

K  *  (L  -  1)  must  be  less  than  15 

Given  the  above  limits  the  rate  1/2,  L  =  7  or  9  code  of  the  modem 
can  be  modeled,  but  the  rate  3/4,  L  =  7  or  9  codes  can  not  be 
modeled.  Another  consideration  for  block  configuration  is  the 
block  requirement  for 

K  *  (samples/bit  in)  =  N  *  ( samples /symbol  out)  (50) 

For  the  rate  1/2  code  with  10  samples /bit  in,  the  coded  bits  have 
only  5  scunples /symbol.  Care  must  be  taken  to  ensure  integer 
relations  exist  or  the  block  will  force  the  result  to  the  nearest 
integer. 

For  timing  purposes,  it  is  noted  that  the  first  valid  output 
from  the  block  is  delayed  according  to  the  equation  [10] 

((K  -  1)  *  (f.,/Rb))  +  (nearest  integer(%T,  *  f„/Rb))  -  1  (51) 
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where 

f„  is  the  sampling  frequency 

Rb  is  the  input  bit  rate 

%T,  is  the  fraction  of  symbol  position:  the  parameter 
that  determines  where  the  input  bit  is  sampled. 

The  clock  output  from  the  block  goes  logic  high  for  the 
first  sample  of  each  output  encoded  bit  after  an  initial  delay 
given  by  Equation  51.  Finally,  a  text  file  must  be  made  to 
specify  the  generator  polynomials  and  values  of  N,  K,  and  L. 
Following  the  required  format  given  in  [10],  with  N=2,  K=l,  L 
=  7,  and  the  modem  generator  polynomials,  the  text  file  code_gen 
was  made. 


Text  File  Code_gen 

2  17 
171 
133 


The  171  in  the  file  is  the  octal  representation  of  the  Go 
generator  polynomial  1111001  (binary).  The  133  is  the  octal 
representation  of  the  Gj  generator  polynomial  1011011  (binary). 
In  addition  to  the  text  file  code_gen,  three  parameters  must  be 
set  in  the  block.  The  seimple  frequency  and  uncoded  bit  rate  are 
required.  Finally,  the  position  in  the  bit  where  the  coder  will 
take  a  sample  is  specified  by  the  fraction  of  symbol  position 
parameter . 
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6.1.3  Viterbi  Decoder  Configuration.  The  SPW™  Viterbi 


decoder  using  soft  decisions  is  well-documented  in  [10].  The 
same  text  file  used  for  configuration  of  the  conventional  encoder 
is  used  by  the  Viterbi  decoder.  The  limits  given  in 
Section  6.1.2  on  N,  K,  L,  also  apply  to  the  decoder.  An 
additional  decoder  limit  is  the  truncation  path  length.  The 
truncation  path  length  is  a  Viterbi  decoder  parameter.  It 
determines  how  much  path  history  will  be  stored  and  used  in  the 
decoder  decision  process.  The  truncation  path  length  limit  is 
100.  The  Viterbi  decoder  implements  soft  decisions  by  first 
quantizing  the  input  data  stream  and  then  making  decisions 
weighted  by  the  confidence  metrics.  A  second  text  file  is 
required  by  the  decoder  block.  This  text  file  provides 
information  on  the  quantization  threshold  levels  and  metrics. 

For  the  modem  model,  eight-level  quantization  of  the  ±1  antipodal 
signal  is  used.  As  determined  from  [6],  the  modem  threshold 
levels  can  be  modeled  as  [  -1,  -0.66,  -0.33,  0,  0.33,  0.66,  1]. 
The  file  named  decode_file  was  created  and  has  the  form  shown 
below. 


Text  File  Decode_file 

8  7  0 

-1.00  6  1 

-0.66  5  2 

-0.33  4  3 

0.00  3  4 

0.33  2  5 

0.66  1  6 

1.00  0  7 
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The  8  in  the  file  specifies  the  number  of  quantization 
levels,  and  the  threshold  levels  are  given  directly  below  the  8. 
The  two  columns  of  numbers  7  through  0  and  0  through  7  are  the 
confidence  levels  given  to  a  received  threshold  level.  The 
second  colvimn  of  numbers,  7  through  0,  is  the  transition  metric 
used  when  a  given  quantization  level  is  received  given  that  a 
logic  low  was  transmitted.  The  third  column  of  numbers,  0 
through  7,  is  the  transition  metric  used  when  a  given 
quantization  level  is  received  given  that  a  logic  high  was 
transmitted.  For  example,  a  quantization  level  1  is  lightly 
weighted  as  a  logic  low  by  the  0  in  column  three,  and  heavily 
weighted  as  a  logic  high  by  the  7  in  column  two.  The  output  of 
the  Viterbi  decoder  is  a  decision  that  a  logic  high  or  low  was 
received  based  on  the  quantization  and  confidence  levels  applied 
to  the  block  input.  In  addition  to  the  two  text  files  code_gen 
and  decode_file,  two  parameters  must  be  set.  The  sample 
frequency  and  encoded  symbol  rate  must  both  be  specified  in  the 
block. 

For  correct  decoding,  the  block  timing  input  must  be  logic 
high  at  the  position  in  the  encoded  bit  where  the  quantization  is 
to  occur  and  logic  low  at  other  times.  The  clock  should  be 
periodic  at  the  encoded  bit  rate.  The  convolutional  encoder 
block  clock  output  is  available  for  timing  the  decoder.  Use  of 
the  encoder  output  clock  for  decoder  timing  makes  the  assumption 
of  perfect  symbol  synchronization  in  the  demodulation  process. 

The  clock  from  the  encoder  must  be  delayed  by  an  amount  equal  to 
any  system  delay  between  the  encoder  output  and  the  decoder 
input.  With  the  encoded  data  input  and  the  correct  timing,  the 


decoder  will  produce  the  original  uncoded  bit  streams  at  the 
block  output.  There  is  a  finite  delay  between  the  first  encoded 
bit  input  to  the  decoder  and  the  first  uncoded  output.  The 
equation  given  in  [10]  for  the  delay  is 

(PL  *  (f„/R.)  *  N)  +  (f3./R.)  *  (N  -  1)  (52) 

where 

f,3  is  the  sampling  frequency 

R,  is  the  input  coded  symbol  rate 

N  is  the  number  of  coded  bits  produced  for  each  uncoded 
bit  input 

PL  is  the  path  length  used  for  Viterbi  decoding 

This  delay  becomes  part  of  the  system  delay  that  must  be 
accounted  for  when  the  original  random  data  at  the  system  input 
is  used  as  a  reference  to  check  the  system  output. 

6.1.4  Back-to-Back  Testing.  The  convolutional  encoder  with 
a  random  data  input  was  connected  directly  to  the  Viterbi  decoder 
in  a  noiseless  environment.  The  encoder  clock  out  was  sent 
directly  to  the  decoder  clock  in,  and  the  inverted  clock  out  was 
sent  to  the  decoder  hold.  The  decoder  hold  will  suspend  decoder 
operation  between  valid  points  in  the  bit  stream.  Signal  sinks 
were  used  to  record  the  random  data  input,  encoded  data,  clock, 
and  decoded  data.  A  test  of  the  system  at  an  uncoded  bit  rate  of 
2400  bps  and  sample  frequency  of  24000  Hz  produced  no  bit  errors. 
The  timing  signal  was  verified  by  inspection  of  the  recorded  data 
as  occurring  at  the  first  sample  of  the  encoded  bit. 
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6 . 2  Coi?plete  System  Testing 


After  the  various  component  blocks  had  been  tested  for 
operation,  the  complete  model  was  assembled.  Figure  14  is  a 
simplified  block  diagram  of  the  system.  The  actual  system  with 


Figure  14.  Complete  System 


signal  sinks  is  shown  in  Figure  38  in  Appendix  A.  The  signal 
sinks  stored  values  that  were  displayed  with  the  Signal  Display 
Editor  (SDE)  part  of  SPW™.  A  constant  one  replaced  the  noise 
input  allowing  verification  of  timing  without  noise.  This 
allowed  system  delays  and  timing  to  be  adjusted  and  verified. 
For  Pb  measurements,  the  signal  sinks  were  removed  and  the  only 
system  output  was  a  file  containing  data  from  the  bit  error 
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counter  block.  The  following  sections  describe  the  system  signal 
flow,  parcuneter  configuration,  timing  considerations,  testing 
process,  and  test  results. 

6.2.1  Signal  Flow.  The  random  data  block  produces  random 
bits  0  for  logic  low  and  1  for  logic  high  at  the  system  bit  rate. 
The  differential  encoder  block  differentially  encodes  the  random 
data,  and  the  output  goes  to  the  convolutional  encoder  block. 

The  convolutional  encoder  generates  a  rate  1/2  constraint  length 
7  encoded  bit  stream.  The  encoded  bit  rate  is  twice  the  uncoded 
bit  rate,  and  the  samples  per  encoded  bit  are  half  the  number  of 
Scunples  per  uncoded  bit.  The  BPSK  modulation  block  produces  a 
baseband  complex  signal  because  the  carrier  parameter  was  set  to 
0  Hz  internal  to  the  block.  The  0  to  1  input  bit  stream  is 
changed  to  an  antipodal  ±1  bit  stream  internally  in  the  BPSK 
modulation  block.  The  modulated  signal  goes  to  a  complex  adder 
and  the  white  noise  given  signal-to-noise  ratio  (SNR)  and 
bandwidth  block.  The  output  of  the  noise  block  is  the  second 
input  to  the  adder.  The  adder  output  is  the  noisy  BPSK  signal 
input  to  the  Costas  loop  demodulator  block.  The  quiescent 
frequency  of  the  voltage  controlled  oscillator  (VCO)  was 
internally  set  to  0  Hz  to  match  the  0  Hz  carrier.  The  Costas 
loop  real_arm  or  in-phase  signal  is  the  demodulated  noisy  signal. 
The  Viterbi  decoder  block  quantizes  the  noisy  input  signal  and 
soft  decision  decodes  the  convolution  encoded  data.  The  encoded 
bit  rate  is  twice  the  decoded  bit  rate,  and  the  number  of  samples 
per  encoded  bit  is  half  the  number  of  samples  per  uncoded  bit  for 
the  rate  1/2  code.  The  differential  decoder  block  differentially 
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decodes  the  Viterbi  decoder  block  output.  The  differential 
decoder  block  output  is  then  input  to  the  bit  error  rate  counter 
block.  The  second  inpuc  to  the  bit  error  rate  counter  block  is 
the  original  bit  stream  after  a  delay  equal  to  the  system  delay. 

6.2.2  System  Parameters .  The  parameter  Sfreq  is  the 
sampling  frequency  and  is  passed  to  all  blocks  except  the  delay 
blocks,  the  inverter  block,  and  the  differential  encode  and 
decode  blocks.  The  SNR  parameter  is  passed  to  the  white  noise 
generator  block.  The  encoded  rate  parameter  is  passed  to  the 
white  noise  block  and  the  Viterbi  decoder  block.  The  parameter 
bit  rate  is  sent  to  the  random  data  block  and  the  convolutional 
encoder  block.  The  loop  frequency,  VCO  gain,  and  arm  filter 
bandwidth  are  all  passed  to  the  Costas  loop.  The  delay  parameter 
is  sent  to  both  the  differential  encoder  and  decoder  blocks.  The 
block  and  initial  error  parameters  are  passed  to  the  bit  error 
rate  counter.  The  T_on  parameter  is  sent  to  both  the  T_on  and 
position  in  symbol  blocks.  Table  13  summarizes  the  parameters 
and  the  parameter  values. 

6.2.3  Timing  Considerations .  The  use  of  the  convolutional 
encoder  clock  output  for  the  Viterbi  decoder  clock  input 
eliminates  symbol  synchronization  concerns  from  the  model.  The 
remaining  timing  considerations  deal  with  timing  the  simulation 
blocks  for  desired  operation.  The  convolutional  encoder  produces 
its  first  valid  output  after  a  delay  determined  by  Equation  51. 
The  uni^  step  block  [l(t-T)]  was  used  to  hold  the  operation  of 
the  BPSK  modulation,  white  noise,  complex  adder,  and  Costas  loop 
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Table  13.  System  Parameters  For  5  Samples /Encoded  Bit  Test 


Parameter 

Name 

Value 

Scimple  frequency 

Sfreq 

24000  Hz 

bit  rate  uncoded 

bitrate 

2400  bps 

bit  rate  coded 

4800  bps 

signal-to-noise  ratio 

SNR 

3.5  dB 

loop  frequency 

loopfreq 

6.5087382  Hz 

VCO  gain 

VCO_gain 

130.17476  Hz/V 

arm  filter  bandwidth 

9600  Hz 

number  of  samples /uncoded  bit 

delay 

-10 

number  of  bits  to  count 

block 

100,000 

first  error  event 

initial_error 

1 

time  to  start  error  count 

T_on 

0.0205  sec 

blocks  until  after  valid  data  was  available  from  the  encoder. 
Another  timing  consideration  was  where  in  the  bit  period  the 
Viterbi  decoder  should  make  a  quantization  sample.  The 
approximate  middle  of  the  bit  was  chosen  to  allow  for  the 
Butterworth  arm  filters  to  settle.  A  bulk  delay  block  was  used 
to  delay  the  clock  output  from  the  convolutional  encoder  by  the 
necessary  eimount  to  produce  a  clock  input  signal  that  caused  the 
Viterbi  decoder  to  sample  the  middle  of  the  input  bits. 

The  final  timing  issues  are  related  to  the  bit  error 
counting  block.  For  the  comparison  of  the  original  data  bits  to 
the  decoded  received  data  bits  to  occur  with  proper  timing,  the 
original  data  bits  must  be  delayed  by  an  amount  equal  to  the 
total  system  delay.  The  Viterbi  decoder  block  introduces  the 
majority  of  the  system  delay.  This  delay  is  determined  by 

I 

!  Equation  52.  Additional  delay  is  produced  by  the  convolutional 

i 
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encoder  as  previously  discussed.  A  bulk  delay  block  is  used  to 
delay  the  original  data  by  the  system  delay.  The  data  bits  input 
to  the  error  counter  are  either  0  for  a  logic  low  or  1  for  a 
logic  high  and  are  ideal  signals.  The  sample  point  for  comparing 
the  bits  is  arbitrary.  The  position  in  symbol  block  by  default 
uses  the  midpoint  of  the  bit  based  on  the  bit  rate  and  sample 
frequency.  The  position  in  symbol  block  is  used  to  allow  the  bit 
error  counter  to  compare  at  the  midpoint  of  the  bit  and  held  the 
block  operation  between  midpoints.  The  T_ON  block  is  used  to 
hold  the  error  counter  block  operation  until  after  the  total 
delay  time.  For  the  system  simulation  with  the  parameter  values 
given  in  Table  13,  Equations  51  and  52  predict  delays  of  4 
Scimples  and  485  samples  respectively.  An  additional  2  sample 
delay  is  caused  by  the  delay  of  the  convolutional  encoder  clock 
output  before  its  use  as  the  clock  input  for  the  Viterbi  decoder. 
The  total  system  delay  is  491  samples  measured  from  the  first 
Scunple  of  the  first  bit  output  from  the  random  data  block  to  the 
first  sample  of  the  first  bit  output  from  the  differential 
decoder  block.  The  value  for  the  T_ON  parameter  is  found  by 
dividing  the  total  system  delay  plus  one  more  sample  by  the 
saunple  frequency.  The  one  additional  sample  starts  the  error 
counter  block  after  the  first  data  bit  transition.  The  value  for 
the  T_ON  parameter  is  492/24000  =  0.0205  sec.  As  recommended  in 
[10],  a  test  run  of  the  system  was  done  with  a  reasonable  number 
of  samples  (5000)  for  this  test.  Additional  signal  sink  blocks 
were  used  to  store  all  timing  signals,  and  the  Signal  Display 
Editor  was  used  to  verify  correct  timing  in  relation  to  the 
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stored  data  bits.  The  full  system  is  shown  in  Figure  39  in 
Appendix  A. 


6.3  Testing  Process 

The  performance  measure  for  the  system  is  probability  of  bit 
error  Pt  given  a  known  energy  per  bit  to  noise  ratio.  The 
white  noise  block  was  used  to  generate  additive  white  Gaussian 
noise  that  was  zero  mean  with  a  variance  adjusted  to  produce  the 
specified  signal-to-noise  ratio.  At  the  point  in  the  system 
where  the  white  noise  block  is  located,  encoded  bits  at  the 
encoded  rate  are  present.  The  white  noise  block  is  operated  with 
the  assumption  that  the  noise  bandwidth  equals  the  bit  rate  of 
the  data  entering  the  block.  As  discussed  in  [4],  the  energy  per 
bit  before  the  rate  1/2  convolutional  encoder  is  3  dB  more  than 
the  energy  per  bit  after  the  encoder.  However,  the  Costas  loop 
arm  filters  were  set  to  twice  the  received  bit  rate  as 
recommended  in  [15].  This  makes  the  noise  bandwidth  twice  what 
the  white  noise  block  operation  assumes.  The  net  effect  is  that 
the  3  dB  more  Et  is  canceled  by  the  3  dB  more  noise  bandwidth. 
This  made  setting  the  SNR  parameter  for  the  white  noise  block  to 
3.5  dB  equivalent  to  having  a  3.5  dB  Eb-to-N„  ratio  into  the 
Costas  loop  demodulator.  With  the  SNR  set,  the  only  remaining 
parameter  to  set  before  the  system  is  ready  for  simulation  is  the 
block  length.  The  block  length  is  determined  by  how  many  bits 
are  to  be  counted  between  error  counter  writes  to  the  results 
file.  By  setting  the  block  length  parameter  to  the  total  number 
of  bits,  for  example  10®,  one  error  count  value  and  one  bit  count 
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value  are  stored  after  each  group  of  10®  bits  are  compared.  For 
testing,  the  block  length  was  set  equal  to  10  divided  by  the 
order  of  magnitude  of  the  expected  error  rate.  The  expected 
error  rates  for  the  system  were  shown  and  discussed  in 
Section  3.2.  The  requirement  for  processing  lO/Pt  bits  was 
discussed  in  Section  3.3.1  No  other  outputs  were  required  from 
the  system,  so  all  signal  sinks  were  removed  from  the  system. 

A  sample  instant  is  considered  an  iteration  when  the 
simulation  is  run.  For  data  bits  with  10  samples  per  bit,  10 
iterations  must  be  processed  for  each  bit.  Therefore,  to 
simulate  10®  bits,  the  simulation  must  run  for  10®  iterations  if 
no  delays  are  present  in  the  simulation.  The  modem  model  system 
delay  was  491  samples.  Each  SNR  ratio  was  simulated  for  a  number 
of  iterations  equal  to  10  times  the  block  length  plus  500 
iterations.  Table  14  shows  the  four  levels  of  SNR,  block  length, 
and  number  of  iterations  used  for  testing. 

Table  14.  Number  of  Iterations  For  5  Samples /Encoded  Bit  Test 


SNR 

Block  Length 

Number  of  Iterations 

3.5  dB 

10® 

1.0005  X  10® 

4.5  dB 

10® 

1.00005  X  10’ 

5.0  dB 

10® 

1.00005  X  10’ 

5.5  dB 

10’ 

1.000005  X  10® 

BPSK  matched  filter  testing.  Section  4.5,  showed  that  over 
one  and  a  half  hours  were  required  to  simulate  10®  bits  with 
1.0001X10*  iterations  in  the  interactive  mode  of  the  Block 
Diagreun  Editor  (BDE).  Therefore,  another  approach  was  taken  for 
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the  modem  model  testing.  The  BDE  tool  for  code  generation 
produces  a  file  containing  a  C  language  program.  The  compiled  C 
code  can  be  run  on  computers  other  than  the  system  with  the  SPW™ 
software,  provided  the  required  SPW^"  library  files  are  available 
on  the  other  computer  [10].  As  a  test  for  consistency  of 
results,  the  BPSK  matched  filter  system  code  was  generated  and 
run.  The  time  for  the  simulation  was  reduced  from  over  an  hour 
and  a  half  to  approximately  25  minutes  with  the  same  resulting 
error  count.  At  the  time  the  C  code  program  is  started,  the 
number  of  iterations  and  noise  seed  value  can  be  entered  on  the 
command  line.  With  a  UNIX  script  file,  the  ten  noise  seeds  for 
each  value  of  iteration  count  can  be  looped  through. 

6.3.1  Initial  Test  Results .  The  simulation  of  the  modem 
model  was  first  done  with  system  parameters  as  shown  in  Table  13 
and  Table  14.  The  number  of  iterations  per  uncoded  data  bit  was 
ten,  whereas  after  the  convolutional  encoder,  there  were  five 
iterations  per  encoded  bit.  The  results  of  the  simulations  are 
shown  in  Table  27  through  Table  30  in  Appendix  B.  Table  15  shows 
the  data  statistics. 


Table  15.  Statistics  For  5  Samples /Encoded  Bit 


a 

3.5  dB 

3.960  X  10-^ 

6.95111  X  10'® 

2.63650  X  10-' 

4.5  dB 

5.314  X  10-' 

3.24627  X  10-® 

5.69760  X  10-® 

5.0  dB 

1.702  X  lO-' 

5.97733  X  10-*° 

2.44486  X  10-® 

5.5  dB 

5.542  X  10'* 

2.89107  X  10-** 

5.37686  X  10-® 
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6«3.2  Five  Samples /Encoded  Bit  Result.  The  mean  value  of 
the  ten  simulation  runs  was  plotted  for  each  Et/N,,  level. 
Additionally,  vertical  lines  representing  the  value  of  ±o  at  each 
Eb/No  level  were  plotted.  In  Figure  15,  the  resulting  curve  is 
labeled  measured  data  (5).  Also  shown,  for  comparison,  are  the 


Figure  15.  Theoretical  and  Measured  Error  Rate  Plots  for 
5  Samples  per  Encoded  Bit 


theoretical  Pb  curves  for  DEBPSK  as  discussed  in  Section  3.1  and 
DEBPSK  with  coding  discussed  in  Section  3.2.  Figure  15  indicates 
that  the  simulation  required  approximately  1  dB  higher  Eb/N^  than 
the  theoretical  line  to  place  Pt  in  the  10"^  range.  The  code  gain 
for  the  simulation  with  five  samples/encoded  bit  is  shown  in 
Table  16.  The  gain  is  measured  as  the  difference  between  the 
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Eb/No  required  to  produce  a  given  Pt  from  Equation  9  (DEBPSK 
theory)  and  the  Eb/N^  of  the  simulation  that  produced  the  given 
Pb. 

Table  16.  Code  Gain  For  5  Samples /Encoded  Bit  Tesc 


Eb/No 

Pb 

Code  Gain 

3.5  dB 

3.960  X  10-^ 

2.6797  dB 

4.5  dB 

5.314  X  lO-”* 

3.2822  dB 

5.0  dB 

1.702  X  10-" 

3.4926  dB 

5.5  dB 

5.542  X  10*® 

3.5993  dB 

6.3.3  Additional  Testing.  The  first  simulations  of  the  full 
system  were  run  with  ten  samples  per  uncoded  bit,  resulting  in 
five  samples  per  coded  bit.  The  effects  of  encoding  and  whether 
the  number  of  saimples  per  bit  in  all  simulation  blocks  should  be 
maintained  at  the  8  to  16  level  were  not  discussed  in  [2].  In 
order  to  determine  if  having  five  samples  per  encoded  bit  was 
affecting  the  system  Pb  performance,  the  system  would  require 
modification  to  achieve  10  samples  per  encoded  bit.  The 
simulation  could  be  modified  in  two  ways  to  achieve  10  samples 
per  encoded  bit.  The  following  sections  discuss  both  methods  of 
modification  and  the  results  of  simulation. 

6. 3 . 3 . 1  Increased  Sample  Frequency.  Increasing  the 
seunple  frequency  in  the  simulation  from  24000  samples/sec  to 
48000  samples/sec  produces  the  desired  10  samples  per  encoded 
bit.  The  system  timing  was  adjusted  to  account  for  the  change  in 
sample  frequency.  The  clock  output  from  the  convolutional 
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encoder  must  now  be  delayed  by  5  samples  to  cause  the  Viterbi 
decoder  to  operate  on  the  midpoint  of  the  encoded  bit.  The  unit 
step  block  was  set  to  transition  on  the  9th  Scimple  as  predicted 
by  Equation  51.  The  total  system  delay  was  calculated  as  in 
Section  6.2.3  to  be  984  samples.  The  parameter  T_ON  remains 
0.0205  sec.  Finally,  the  delay  parameter  passed  to  the 
differential  encoder  and  decoder  was  changed  to  20  because  there 
were  20  samples  per  uncoded  bit.  Figure  40  in  Appendix  A  sh-.ws 
the  system.  The  C  code  for  the  modified  system  was  generated  and 
the  simulation  run  for  a  SNR  of  3.5  dB.  Table  32  in  Appendix  B 
shows  the  results.  Table  17  shows  the  statistics  of  the  results. 

Table  17.  Statistics  For  48000  Hz  Sample  Frequency; 

3.5  dB  Efa/No 


n  (10'^) 

(10'®) 

a  (10'®) 

1.924 

4.62933 

2.15159 

Although  simulation  time  was  not  a  area  of  study  for  this  thesis, 
it  became  a  concern.  With  20  samples  per  uncoded  bit,  the  number 
of  iterations  for  the  SNR  of  3.5  dB  simulation  was  2.001X10®. 

Over  one  and  a  half  hours  were  required  for  each  noise  seed 
simulation  if  the  process  had  sole  use  of  the  computer. 

6 . 3 . 3 . 2  Decrease  Bit  Rate.  Instead  of  raising  the 
saunple  frequency,  the  data  rate  can  be  lowered  to  reach  10 
seunples  per  encoded  bit.  The  bit  rate  parameter  change  from 
2400  bps  to  1200  bps  caused  the  need  to  change  the  arm_band 
pareuneter  from  9600  Hz  to  4800  Hz.  System  delays  were  changed  to 
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the  same  values  that  resulted  from  doubling  the  sample  frequency 
as  discussed  previously.  The  T_ON  parameter  was  changed  to 
0.041  sec  to  account  for  the  fact  that  at  1200  bps  it  took  twice 
as  long  to  produce  the  same  number  of  bits  as  were  produced  at 
2400  bps.  The  modified  system  is  shown  in  Figure  41  in 
Appendix  A. 

The  C  code  for  the  modified  system  was  generated  and  the 
simulation  run  for  a  SNR  of  3.5  dB.  Again,  2.001X10®  iterations 
were  processed  to  produce  1X10®  data  bits.  Table  33  in 
Appendix  B  shows  the  results.  Table  18  shows  the  statistics  of 
the  results. 

Table  18.  Statistics  For  1200  bps  Data  Rate;  3.5  dB  Et/No 


H  (10-^) 

(10*®) 

o  (10*®) 

1.842 

8.33289 

2.88667 

6 . 3 . 3 . 3  Ten  Samples /Encoded  Bit  Test.  The  two  methods 
of  simulation  with  10  samples  per  encoded  bit  through  the  Costas 
loop  produced  nearly  the  same  probability  of  bit  error  for  a  Eb/N^ 
of  3.5  dB.  The  difference  between  the  probability  of  bit  error 
result  for  each  method  was  less  than  the  value  of  a  for  either 
method.  Based  on  this  fact,  the  remaining  simulations  were  run 
with  the  system  that  achieved  20  samples  per  uncoded  bit  by 
lowering  the  data  rate.  An  added  consequence  of  using  this 
method  is  that  the  Costas  loop  is  configured  as  it  was  for 
testing  described  in  Section  5.5. 
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The  simulation  was  run  for  each  of  the  10  noise  seeds  at  the 


four  SNR  levels.  Table  19  shows  the  four  SNR  levels,  block 
length,  and  number  of  iterations.  Table  20  shows  the  system 
parcimeters . 


Table  19.  Number  of  Iterations  for  10  Samples /Encoded  Bit  Test 


SNR 

Block  Length 

Number  of  Iterations 

3.5  dB 

10^ 

2.001  X  10® 

4.5  dB 

10® 

2.0001  X  10’ 

5.0  dB 

10® 

2.0001  X  10’ 

5.5  dB 

10’ 

2.00001  X  10® 

Table  20.  System  Parameters  For  10  Samples/Encoded  Bit  Test 


Parameter 

Name 

Value 

sample  frequency 

Sfreq 

24000  Hz 

bit  rate  uncoded 

bitrate 

1200  bps 

bit  rate  coded 

2400  bps 

signal-to-noise  ratio 

SNR 

3.5  dB 

loop  frequency 

loopfreq 

6.5087382  Hz 

VCO  gain 

VC0_gain 

130.17476  Hz/V 

arm  filter  bandwidth 

arm_band 

4800  Hz 

number  of  samples/uncoded  bit 

delay 

-20 

number  of  bits  to  count 

block 

100,000 

first  error  event 

1 

time  to  start  error  count 

T_on 

0.041  sec 
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The  simulation  time  using  compiled  C  code  ranged  from 
approximately  1.5  hours  for  2.001X10®  iterations  to  about  27 
hours  for  2.0001X10®  iterations.  The  times  given  above  are  for  a 
single  noise  seed  and  a  SPARC-2”*  computer.  With  three  computers 
running  simulations  with  different  noise  seeds,  over  four  days  of 
computer  time  were  required  to  simulate  the  system  at  the  5.5  dB 
SNR  level  for  all  10  noise  seeds.  The  simulation  results  are 
shown  in  Table  35  through  Table  38  in  Appendix  B.  Table  21  shows 
the  statistics  of  the  results. 

Table  21.  Statistics  For  10  Samples /Encoded  Bit 


Efa/N, 

a 

3.5  dB 

1.842  X  10'^ 

8.33289  X  10-® 

2.88667  X  10-® 

4.5  dB 

2.332  X  10-* 

6.75733  X  10-*° 

2.59949  X  10-® 

5.0  dB 

7.220  X  10*® 

2.47511  X  10-*° 

1.57325  X  10-® 

5.5  dB 

1.986  X  10-® 

4.49822  X  10-*^ 

2.12090  X  10-* 

6.3.4  Ten  Samples /Encoded  Bit  Result.  The  mean  value  of  the 
ten  simulation  runs  was  plotted  for  each  Et/N^  level.  Vertical 
lines  representing  the  value  of  ±a  at  each  E^/No  level  were  also 
plotted.  In  Figure  16,  the  resulting  curve  is  labeled  measured 
data  (10).  The  other  curves  are  the  same  as  shown  in  Figure  15. 
The  simulation  with  10  samples  per  encoded  bit  requires 
approximately  0.7  dB  higher  Eb/N^  than  the  theoretical  line  to 
place  Pb  in  the  10*®  range.  Table  22  gives  the  code  gain  for  this 
simulation. 
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Table  22.  Code  Gain  For  10  Samples /Encoded  Bit  Test 


Eb/No 

Pb 

Code  Gain 

3.5  dB 

1.842  X  10*^ 

3.3570  dB 

4.5  dB 

2.332  X  10-^ 

3.8068  dB 

5.0  dB 

7.220  X  10*® 

3.9633  dB 

5.5  dB 

1.986  X  10*® 

4.0914  dB 

6.4  Results  Discussion 

The  improvement  in  the  simulation  performance  can  be  seen  in 
Figure  16.  The  10  samples  per  encoded  bit  is  approximately 
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0.4  dB  closer  to  the  theoretical  line  than  the  result  of  testing 
with  5  samples  per  encoded  bit.  The  theoretical  signal-to-noise 
ratio  for  the  loop  (p^)  ranged  from  approximately  19.5  d3  to 

21.5  dB  as  discussed  in  Chapter  V.  As  shown  by  [5],  and 
discussed  in  Chapter  III,  the  Pb  vs  Et/N^  curve  for  this  pi  should 
be  steeper  and  differ  from  the  theoretical  case  of  p^  =  oo  by  less 
than  0.5  dB.  The  Pb  curve  resulting  from  the  simulation  appears 
to  agree  with  a  p^  on  the  order  of  13.5  dB.  The  actual  pi  of  the 
loop  simulation  was  not  found  by  the  testing  done  in  this  thesis. 
Simulations  that  would  allow  determination  of  the  actual  Costas 
loop  Pi  must  be  performed  before  an  exact  determination  of  the 
effect  of  Pi  on  the  results  would  be  possible. 

6.5  Summary 

This  chapter  first  covered  the  configuration  and  testing  of 
the  coding  blocks.  The  ability  of  the  coding  blocks  to  encode 
and  decode  a  bit  stream  without  error  in  a  noiseless  test 
verified  the  correct  configuration  of  these  blocks.  The  full 
model  was  then  described  with  a  detailed  discussion  of  timing 
considerations.  The  majority  of  the  system  delay,  over  98%,  was 
the  result  of  the  convolutional  encoder  and  Viterbi  decoder 
blocks.  The  initial  test  result  was  improved  by  approximately 
0.4  dB  by  doubling  the  number  of  samples  per  encoded  bit 
processed  by  the  Costas  loop.  Both  increasing  the  sample 
frequency  or  decreasing  the  data  rate  doubled  the  number  of 
samples  per  encoded  bit  and  resulted  in  an  improvement  in  Eb/No 
required  to  achieve  a  given  Pb.  The  full  modem  simulation 
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resulted  in  a  code  gain  of  between  3.357  dB  and  4.0914  dB  for 
Efa/No  values  from  3.5  dB  to  5.5  dB.  These  results  differ  from  the 
theoretical  code  gain  by  approximately  0.7  dB  for  Pt  in  the  10** 
range.  It  is  unknown  if  pi  contributed  to  the  difference  from 
theoretical,  because  testing  performed  did  not  calculate  the 
actual  Pi  of  the  simulated  Costas  loop.  Before  these  results  are 
used  as  more  than  a  rough  estimate  of  the  CQM-248  performance, 
further  refinement  of  the  simulation  should  be  done.  Additional 
testing  of  the  Costas  loop  and  refinements  to  the  simulation  are 
proposed  in  Section  7.3. 
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VII.  Conclusions  and  Recommendations 

7. 1  Conclusions 

The  goal  of  this  thesis  was  to  estimate  the  CQM-248  modem 
performance  by  computer  simulation.  The  model  of  the  modem 
includes  differential  and  convolutional  encoding  of  the  data  with 
BPSK  modulation  of  the  encoded  data.  The  channel  was  modeled  as 
a  AWGN  baseband  disturbance.  A  second-order  Costas  loop 
demodulates  the  received  noisy  data  and  soft  decision  Viterbi 
decoding  is  used.  A  differential  decoder  provides  the  final  data 
output  for  comparison  to  the  original  data.  Accomplishing  this 
goal  required  three  objectives. 

The  first  objective  was  to  gain  experience  with  the  SPW™ 
tool.  The  BPSK  matched  filter  detection  system  was  used  to  gain 
experience  with  the  SPW™  tool.  Also,  testing  the  matched  filter 
system  demonstrated  the  Monte  Carlo  error  estimation  method. 

Most  importantly,  the  matched  filter  system  verified  the 
operation  of  the  bit  error  counter  and  white  noise  generation 
blocks . 

The  matched  filter  detection  system  proved  useful  as  a 
learning  tool.  Simulation  timing  issues  were  discovered  and 
solved.  For  example,  the  delay  of  the  original  data  for 
comparison  with  the  received  data  was  not  obvious.  The  system 
worked  as  expected  with  P^  estimated  results  within  one  standard 
deviation  of  theoretical. 

The  second  objective  was  to  develop  and  demonstrate  proper 
performance  of  a  Costas  loop  demodulator.  The  loop  response  to  a 
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phase  step  or  a  frequency  step  was  simulated  to  verify  loop 
operation  in  a  noiseless  case.  The  loop  response  was  compared  to 
theoretical  predictions. 

The  Costas  loop  with  the  corrected  equations  in  the  loop 
filter  was  used  in  further  simulations.  The  loop  rise  time  to 
peak  deviation  in  response  to  a  phase  step  was  within  2.9%  of 
theoretical.  The  peak  deviation  in  response  to  a  phase  step  was 
within  7.7%  of  theoretical.  For  a  step  in  frequency,  the  loop 
steady  state  phase  error  was  within  5.5%  of  theoretical,  and 
lock-in  occurred  within  one  cycle  of  the  loop  natural  frequency. 
This  comparison  to  theoretical  predictions  demonstrated  loop 
performance  in  a  noiseless  test. 

The  final  objective  was  to  assemble  and  test  a  full  system. 
The  complete  model  was  assembled  after  the  Costas  loop  was 
correctly  working.  The  full  system  included  differential 
encoding  and  decoding;  convolutional  encoding  with  eight-level 
soft  Viterbi  decoding;  BPSK  modulation  of  the  encoded  data;  an 
AWGN  channel;  and  the  Costas  loop  demodulator.  Two  ratios  of 
samples  per  bit  were  used  in  the  testing  and  four  Eb/N^  levels 
were  tested. 

The  full  modem  simulation  resulting  Pb  was  improved  by 
approximately  0.4  dB  by  doubling  the  number  of  samples  per 
encoded  bit  processed  by  the  Costas  loop  from  5  samples  per 
encoded  bit  to  10  samples  per  encoded  bit.  However,  the  Pb  curve 
still  differed  from  the  theoretical  curve  by  about  0.7  dB  for  Pb 
in  the  10'^  range.  As  discussed  in  Chapter  III  and  Chapter  VI, 
the  loop  signal  to  noise  ratio  affects  the  slope  of  the  Pb  curve. 
Because,  the  actual  loop  signal-to-noise  ratio  was  not  determined 
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for  the  simulation,  the  effect  on  the  simulation  result  is 
unknovm.  The  full  modem  simulation  resulted  in  a  code  gain  of 
between  3.36  dB  and  4.09  dB  for  Eb/N,,  values  from  3.5  dB  to 
5.5  dB.  Further  refinement  of  the  simulation  should  be  done 
before  these  results  are  used  as  more  than  a  rough  estimate  of 
the  CQM-248  performance. 

7.2  Recommendations 

The  following  paragraphs  are  recommendations  for  follow-on 

work. 

1.  Testing  of  the  Costas  loop  should  be  done  in  a  noisy 
environment.  The  simulated  Costas  loop  signal  to  noise  ratio 
should  be  determined  and  compared  to  the  theoretical  predictions 
in  Chapter  V.  This  comparison  could  lead  to  further 
modifications  to  the  Costas  loop  that  would  affect  Pt  of  the  full 
system  simulation. 

2.  The  model  should  be  extended  to  a  coded  QPSK  system. 
This  would  provide  another  measure  of  modem  performance.  The 
building  block  for  the  extension  to  QPSK  would  be  the 
modification  of  the  Costas  loop  to  simulate  a  four-phase  Costas 
loop.  Phase  and  frequency  step  responses  of  the  Costas  loop 
should  again  be  simulated  and  compared  to  theoretical 
predictions . 

3.  When  the  CQM-248  modem  is  available,  hardware  testing 
should  be  done.  The  Pj,  vs  Eb/N^  results  should  then  be  compared 
to  the  simulation. 


88 


The  plots  included  in  this  appendix  were  generated  by  the 
Block  Diagrcun  and  Signal  Display  Editor.  Figures  17  through  25 
and  Figures  34  through  41  were  generated  by  the  Block  Diagram 
Editor  tool  from  SPW™.  Figures  26  through  33  were  created  by  the 
Signal  Display  Editor  tool  from  SPW™‘  Figure  17  is  the  BPSK 
matched  filter  system.  Figures  18-20  are  the  detail  block 
diagrcuns  used  in  the  matched  filter  system.  Figures  21-25  are 
the  Costas  loop  detail  and  the  diagrams  of  the  systems  used  to 
test  the  loop.  Figures  26-33  are  the  Costas  loop  test  system 
outputs.  Figures  34  and  35  are  the  detail  and  symbol  for  the 
differential  encoder.  Figures  36  and  37  are  the  detail  and 
symbol  for  the  differential  decoder.  Figures  38-41  are  block 
diagrams  of  the  systems  used  for  full  system  testing. 
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Figure  17.  Block  Diagram  of  BPSK  Matched  Filter  System 


Figure  18.  White  Noise  Generator  Detail 


MflTCIIED.riLTER  UEMOUULniOR  BLOCK  PYROMETERS 


Figure  19.  PSK  Matched  Filter  Demodulator  Detail 


PSK  DETECTOR  BLOCK  PPRmETERS 


Figure  20.  PSK  Detector  Detail 


Modified  Costas  Loop  Detail 


SB  1  500 
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23. 


21000.0  Hz 


Figure  24.  Costas  Loop  Phase  Step  System 


Figure  25.  Costas  Loop  Frequency  Step  Test  System 


to  Phase  Step  Plot 


Figure  29.  Costas  Loop  Response  to  Phase  Step  Plot 
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Figure  31.  Costas  Loop  Response  to  Frequency  Step  Plot 


Figure  32.  Costas  Loop  Response  to  Frequency  Step  Plot 


Figure  33.  Costas  Loop  Response  to  Frequency  Step  Plot 


Figure  34.  Differential  Encoder  Detail 
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Figure  35.  Differential  Encoder  Symbol 


Differential  Decoder  Detail 
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Figure  37.  Differential  Decoder  Symbol 


Ill 


Figure  38.  Full  System  for  Timing  Verification 
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Figure  39.  Full  System  Used  for  Five  Samples  per  Coded  Bit  Test 
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at  1200  Bit  per  Second 


Appendix  B 

Table  23  through  Table  38  are  the  data  results  and 
statistics  for  estimation  of  probability  of  bit  error  as  found 
using  the  Monte  Carlo  Method. 

Table  23.  Matched  Filter  Detection  BPSK  For  4  dB  Eb/N^ 

Run  Pb  (10~^) 

_ 1 _ 2.10 _ 

_ 2 _ 1.20 _ 

_ 3 _ 1.90 _ 

_ 4 _ 1.40 _ 

_ 5 _ 1.00 _ 

_ 6 _ 1.10 _ 

_ 7 _ 0.60 _ 

_ 8 _ 1.20 _ 

_ 9 _ 0.90 _ 

10  2.00 


Table  24.  Pb  Matched  Filter  Detection  BPSK  For  6  dB  Eb/No 

Run  Pb  (10*^) 


2. 

40 

3. 

60 

3. 

20 

3. 

50 

3. 

20 

1. 

70 

3. 

20 

2. 

50 

2. 

90 

2. 

20 

116 


(10~M 

1.842 


(10-^) 

8.33289 


L 


o  (10~*) 
2.88667 
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